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Abstract 

Precision predictions for phenomenologically interesting observables such as the tt in- 
variant mass distribution and forward-backward asymmetry in top-quark pair production 
at hadron colhders require control over the differential cross section in perturbative QCD. 
In this paper we improve existing calculations of the doubly differential cross section in 
the invariant mass and scattering angle by using techniques from soft-collinear effective 
theory to perform an NNLL resummation of threshold logarithms, which become large 
when the invariant mass M of the top-quark pair approaches the partonic center-of-mass 
energy We also derive an approximate formula for the differential cross section at 
NNLO in fixed-order perturbation theory, which completely determines the coefficients 
multiplying the singular plus distributions in the variable (1— M^/s). We then match our 
results in the threshold region with the exact results at NLO in fixed-order perturbation 
theory, and perform a numerical analysis of the invariant mass distribution, the total 
cross section, and the forward-backward asymmetry. We argue that these are the most 
accurate predictions available for these observables at present. Using MSTW2008NNLO 
parton distribution functions (PDFs) along with as{Mz) = 0.117 and mt = 173.1 GeV, 
we obtain for the inclusive production cross sections at the Tevatron and LHC the values 
CTcvatron = (6-30 it O.lO^'^Qgs) pb and cjlhc = (149±7±8) pb, where the first error results 
from scale variations while the second reflects PDF uncertainties. 



1 Introduction 



The top quark is the heaviest known particle in the Standard Model (SM) of fundamental 
interactions. Because of its large mass, it is expected to couple strongly with the fields re- 
sponsible for electroweak symmetry breaking, and the detailed study of top-quark properties 
is likely to play a key role in elucidating the origin of particle masses. The measurement 
of top-quark related observables is therefore one of the main goals of the Fermilab Tevatron 
and CERN Large Hadron Collider (LHC). To date, thousands of top-quark events have been 
observed by two different experiments at the Tevatron, and the top-quark mass has been ex- 
tracted at the percent level [1] . The experiments at the LHC are expected to observe millions 
of top-quark events per year already in the initial low-luminosity phase, bringing the study of 
the top-quark properties into the realm of precision physics. In particular, the total inclusive 
cross section for top-quark pair production is expected to be measured with a relative error of 
5% to 10% at the LHC [2]. In addition to the mass and the total inclusive cross section, kine- 
matic distributions and charge asymmetries are also of interest. For instance, the tt invariant 
mass distribution can be used as a complementary method for measuring rrit [3] . The presence 
of bumps in the smoothly decreasing ti invariant mass distribution would be a clear signal of 
an s-channel heavy resonance, which is predicted in many new physics scenarios [3-5]. Such 
searches have been pursued at the Tevatron [6-8], and results for the top-pair invariant mass 
distribution were recently obtained from data collected by the CDF collaboration [9]. The 
forward-backward asymmetry in tt production has also been measured recently [10-12] and 
shows a potential deviation from the theoretical predictions. 

To make optimal use of the experimental measurements requires theoretical predictions 
of similar precision. At hadron colliders, a leading-order (LO) prediction is usually insuf- 
ficient, and including higher-order corrections in QCD is mandatory. Current theoretical 
predictions for ti production in fixed-order perturbation theory are based on next-to-leading 
order (NLO) calculations of the total cross section [13-16], differential distributions [17-19], 
and the forward-backward asymmetry [20,21]. The NLO computations suffer from theoretical 
uncertainties larger than 10%, both for Tevatron and LHC center-of-mass energies. These 
uncertainties are partly due to our imperfect knowledge of the parton distribution functions 
(PDFs), and partly to the truncation of the perturbative series in the strong coupling constant, 
which introduces a dependence on the unphysical renormalization and factorization scales into 
the physical predictions. The latter uncertainty is typically reduced by including more terms 
in the perturbative series, and for this reason the calculation of the partonic cross section 
to next-to-next-to-leading order (NNLO) has been an area of active research. This includes 
studies of the two- loop virtual corrections [22-26] , the squared one- loop corrections [27-29] , 
and the NLO corrections to the final state ti + jet [30]. However, due to the complexity of 
the calculations, complete results for the cross section at NNLO are not yet available. 

In the absence of full NNLO results, one way to improve the NLO predictions is to utilize 
threshold resummation methods [31,32] to incorporate some (presumably) dominant contri- 
butions from higher orders. For the total cross section, a common approach is to work in 
the limit where the top-quark pair is produced nearly at rest, and to resum logarithms in the 
velocity (3 = y^l — 4m^ / s — )■ 0, with s the partonic center-of-mass energy. Such resumma- 
tions have been carried out at leading- logarithmic (LL) [33-38], next-to-leading- logarithmic 
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(NLL) [39], and approximate next-to- next-to-leading-logarithmic (NNLL) order [40-42]. Only 
very recently have the complete expressions at NNLL order been obtained [43-45]. The calcu- 
lation of differential cross sections is more involved due to the appearance of several kinematic 
variables. A typical approach in that case is to work in the threshold limit where the parame- 
ter (1 — 2;) = 1 — M^/ s — )■ 0, with M the invariant mass of the ti pair, and to resum logarithms 
in (1 — 2;). In that case real gluon emission is soft, but the parameter /3 is a generic 0{1) 
variable. The current frontier for resummation in this limit is NLL calculations for differential 
cross sections [46-48] and the forward-backward asymmetry [49]. In addition, approximate 
NNLL calculations were performed in [50-53]. Extending these results to full NNLL order 
has been made possible by our recent calculation of the two-loop anomalous-dimension ma- 
trices [54,55]. We have presented an approximate NNLO formula for the ti invariant mass 
distribution in [56]. The goal of the present paper is to derive a renormalization-group (RG) 
improved expression for the doubly differential cross section at NNLL order, in which all 
threshold-enhanced terms are resummed. We will match this expression with the exact fixed- 
order NLO results and study the top-pair invariant mass distribution, the forward-backward 
asymmetry, and the total cross section at NLO-I-NNLL order. An important part of our anal- 
ysis is a comparison of the total cross section obtained here with those obtained in the limit 
/3 — )■ 0. The outcome of this comparison leads us to conclude that the predictions obtained in 
this paper are the most precise available at present. 

The paper is organized as follows. In Section 2 we review the kinematics and the struc- 
ture of factorization in the threshold region. We then derive the factorization formula for the 
hard-scattering kernels into products of hard and soft matrices using soft-coUinear effective 
theory (SCET) in Section 3. In Section 4 we present the calculation of the hard and soft 
matrices at NLO, and describe several checks on our results. Section 5 deals with the RG 
properties of the hard and soft functions. We derive a formula for the resummed cross sec- 
tion in momentum space using RG methods and describe its evaluation at NNLL order. We 
also review the derivation of the approximate NNLO formula, which has been presented first 
in [56]. In Section 6 we perform numerical studies of the invariant mass distribution, the total 
cross section, and the forward-backward asymmetry, utilizing both RG-improved perturbation 
theory at NNLL order and the NNLO approximate formula. This section includes the afore- 
mentioned comparison of different methods for obtaining the total cross section. We conclude 
in Section 7. 

2 Kinematics and factorization at threshold 

We consider the process 

NiiPi) + N2{P2) ^ t{p^) + i{p,) + X{px) , (1) 

where X is an inclusive hadronic final state ^. At Born level this proceeds through the qq 
annihilation and gluon-fusion channels 

q{pi) + q{p2) tips) + t{p4) , 

^Throughout the analysis we treat the top-quarks as on-shcU partons and neglect their decay. Corrections 
to this picture are suppressed by Tt/mt ^ 1. 
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9{Pi)+9{P2) ^t{ps)+t{p^), (2) 
where pi = XiPi and p2 = X2P2- We define the kinematic invariants 

S={Pi + P2f , S = {pi+P2f , = {P3+P,y , 

ti = {Pi - Pzf - ml , Ml = {p2 - pzf - ml , (3) 

and momentum conservation at Born level implies s + ti + Mi = 0. 

In this section we consider the structure of the differential cross section near the partonic 
threshold. While the fully differential cross section depends on three kinematic variables, in 
this paper we are mainly interested in the doubly differential cross section expressed in terms 
of the invariant mass M of the ti pair and the scattering angle 9 between pi and p^ in the 
partonic center-of-mass frame. To describe this distribution we introduce the variables 



The quantity [5t gives the 3- velocity of the top quarks in the ti rest frame. One often considers 



a related variable /3 = y 1 — 4ml /s, "which coincides with (3^ in the limit z — )■ 1. 

Before moving on, we would like to clarify our definition of the so-called "threshold region" . 
In general, the definition of this region depends on the measured observable. For instance, 
to describe the total inclusive cross section, one often considers the limit /3 — [39-45]. To 
describe the ti invariant mass distribution, on the other hand, we consider the limit z ^ 1, 
with a generic C(l) variable. In that case the top quarks are not produced at rest, nor are 
they highly boosted, and the emitted partons in the final state X are constrained to be soft. 
This threshold limit is often referred to in the literature as the pair-invariant-mass kinematics, 
and the theoretical framework to deal with this situation was developed in [46,47,57,58]. 
A closely related framework can be applied to processes where only a single top-quark is 
observed [59,60], but will not be considered here. Threshold resummation at fixed invariant 
mass M is completely analogous to the well-studied case of Drell-Yan or Higgs production, 
and it becomes more and more important as M is increased. Therefore, the RG-improved 
result for da/dM offers the best possible prediction for this observable. Moreover, integrating 
this distribution over M then yields the best possible prediction for the total cross section. In 
particular, we believe that in this way we obtain a better description of the total cross section 
than the usual approach of resumming soft gluon and Coulomb terms for /3 — )■ 0. The reason 
is that at the Tevatron and LHC, the cross section is dominated by events with f3 ~ 0.3-0.6, 
which is not especially small. We will come back to this important point in Section 6.3. 

According to the QCD factorization theorem [61], the differential cross section in M and 
cos 6 can be written as 

d'^a 87r/3j ^ dz 



dMd cos 9 3sM 



/" dz 

Yl / — ffiA^/^^l^f)Cij{z,M,mt,cose,iXf) , (5) 



where /i/ is the factorization scale, and the parton luminosity functions ffij are defined by 



"1 dx 

iijiy,f^.f)=l — fi/Nr{x,flf) fj/N^{y/x,Hf) . (6) 



y 



X 
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The /j/AT are universal non-perturbative PDFs of the parton i in the hadron iV, which can be 
extracted from experimental data. The hard-scattering kernels Cij are related to the partonic 
cross sections and can be calculated as a power series in as- We shall write their expansion as 



At leading order in as, only Cqq and Cgg are non-zero. They are proportional to 5(1 — z) 
read 



(7) 



and 
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(8) 



where = 3 is the number of colors in QCD, and ti and ui can be expressed in terms of M 
and cos 6' as 



tl 



(1 — f3t cos 9) , Ul 



M2 



/3f cos 6) . 



(9) 



Note that for the doubly differential cross section the coefficient Cqq is also needed. It can be 
obtained from the expression for Cqq by replacing cos 6* — t- — cos 6, which is a symmetry at 
tree level but not beyond. 

At higher orders in ag the hard-scattering kernels receive corrections from virtual loop 
diagrams and real gluon emissions in the qq and gg channels, as well as from other partonic 
channels such as gq — > tiq. The calculation of these corrections near threshold is greatly 
simplified. For z ^ 1 there is no phase-space available for hard gluon emission, which is thus 
suppressed by powers of {1 — z). Moreover, contributions from channels such as gq — >■ ttq, 
which involve external soft-quark fields, are also suppressed. The partonic scattering process 
is thus dominated by virtual corrections and the real emission of soft gluons. The phase-space 
for such processes is effectively that for a two-body final state, so the hard-scattering kernels 
can be written in terms of the kinematic invariants from the Born-level processes. Therefore, 
up to corrections of order {1 — z), we can rewrite (5) in the threshold region as 



dMd cos e 3sM 



^ dz 

z 



igg\ 



-/ z, fif) Cgg{z, M, nit, COS 9,fif) 



+ 



Iqq 



(r/ z, fif) Cqq{z, M, mt, cos 6, /i/) + #gg(r/ z, fif) Cqq{z, M, rrit, - cos d, fif) , (10) 



where ffqq is understood to be summed over all light quark fiavors. 

In the threshold limit z ^ 1, the hard-scattering kernels Cij can be factorized into a 
product of hard and soft functions according to 



Cij{z, M, rrit, cos6, fif) = Tr Hij{M, rrit, cos6, fif) Sij{\/§{1 — z), rrit, cos6, fij) + 0{1 
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The boldface indicates that the hard functions i/y and soft functions are matrices in color 
space, with respect to which the trace is taken. We will derive this formula in the next section, 
using techniques from SCET. (A similar factorization formula in Mellin moment space was 
derived in [47].) The hard functions are related to the virtual corrections and are ordinary 
functions of their arguments, while the soft functions are related to the real emission of soft 
gluons and contain singular distributions in {1 — z). In addition to terms proportional to 
6{1 — z), the n-th order corrections in as also contain plus distributions of the form 



^"'(l - z) 
I- z 



m 



0, 



2n - 1 



(12) 



where 



dz 



in^ri 



dz^ '-[9{z)-g{l)]-g{l) dz ^ (13) 

1- z Jq 1- z 



for an arbitrary function g{z). These singular distributions make the perturbative series badly 
convergent near threshold and must be resummed to all orders in perturbation theory. In this 
paper we perform such a resummation directly in momentum space [62], up to NNLL order. 
To this end, we extend the procedure for deep-inelastic scattering [63], Drell-Yan process [64], 
Higgs production [65,66], and direct photon production [67] to processes with four colored 
external particles. The formalism will be described in Section 5. 



3 Factorization at threshold in SCET 

In this section we derive the factorization formula (10) for the hard-scattering kernels in the 
threshold region using SCET [68-70] and heavy-quark effective theory (HQET) (for a review, 
see [71]). The derivation is similar to the ones in [63,64,67], but is more complicated due to 
the presence of additional Dirac and color structures. The derivation of factorization in the 
effective theory relies on a two-step matching procedure. In the first step, fluctuations at the 
hard scale from virtual corrections are integrated out by matching QCD onto an effective- 
theory with coUinear and soft degrees of freedom. The Wilson coefficients from this matching 
step give the hard function when squaring the amplitude. In the second step, the soft degrees 
of freedom are integrated out, giving rise to a soft function, which is defined as the vacuum 
expectation value of a Wilson loop operator. 

3.1 Fields and operators 

The scattering amplitude for tt production involves several scales, which we assume to satisfy 

s,M\ > s(l - zy > A^cD (14) 

in the threshold region. The elements of the first set of scales are taken to be of the same order 
and shall be collectively referred to as hard scales, whereas s{l — z)'^ defines the soft scale. The 
small quantity \ = {1 — z) <til then serves as the expansion parameter in the effective theory. 
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Note that we treat M and rut as of the same order, which means that the top quarks are not 
highly boosted. To describe the invariant mass spectrum in the region where M ^ 2mt, a 
more appropriate treatment would require a different effective theory to separate these two 
scales, and two jet functions have to be introduced for the top and anti-top quarks. Such an 
approach was adopted in [72] for top-quark production in e~^e~ collisions, where rrit/M was 
used as a small expansion parameter. However, given that even for M as large as 1.5 TeV the 
ratio 2mt/M ^ 0.23 is still a reasonable 0{1) parameter, we see no need to adopt it for the 
present work. 

The formalism for SCET applied to a generic ra-body scattering process involving both 
heavy and light partons was set up in [73,74]. In our case, the effective theory contains two 
sets of collinear fields to describe the degrees of freedom in the incoming hadrons, two sets of 
HQET fields to describe the outgoing heavy quarks, and a single set of soft fields describing 
the final state X and the soft interactions among particles. In classifying the collinear fields 
we define two light-like vectors n and n in the directions of the colliding partons, which satisfy 
n ■ n = 2. The collinear quark fields are related to the QCD fields by 

^^{x) = ^i;{x), e^(a;) = ^^(x). (15) 

The collinear gluon fields in a single collinear sector are identical to those in QCD, with their 
momenta restricted to be collinear to the given direction. In constructing operators below, it 
will be convenient to introduce the manifestly gauge-invariant combinations of fields [69,75] 

Xn{x) = Wlix) e„(x) , <j_(x) = Wlix) [iD^^Wnix)] , (16) 

where the n-coUinear Wilson line is defined by 

Wn{x) = V exp J ds n ■ An{x + sn)^ , (17) 

and V denotes path ordering. The corresponding objects for the ra-coUinear fields are obtained 
by interchanging n and n. The HQET fields hy,^ and h^j^ are labeled by the velocities of the 
top quark and anti-quark, which are related to their momenta as 

Ps = mtvi^ + ^3 , P4= ^tVi + K ■ (18) 

The residual momenta ki scale as soft momenta and are set to zero for on-shell quarks. In 
terms of the QCD top-quark fields, the HQET fields are defined as 

KX^) = '^e-'^'''-''t{x). (19) 

A crucial property of the leading-order SCET Lagrangian is that the interactions of soft 
gluon fields with collinear and heavy-quark fields are described by eikonal vertices. The explicit 
form of the interaction terms for soft gluons with the fermion fields is 

Ant = in{x)^gn ■ As{x)in{.x) +ln{.x)^gn- As{x)in{.x) 
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+ hy^ (x) gv3-A, (x) hy^ (x) + V (x) gv^- (x) V i^) > (20) 

and those between collinear and soft gluon fields can be deduced by making tfie substitution 
An An + n ■ Asn/2 (and similarly for the n-coUinear fields) in the Yang-Mills Lagrangian. 
Such eikonal interactions can be absorbed into Wilson lines via the field redefinitions [69, 76] 

x;:(x) ^ [Sn{x)Y' xT\x) , ^J:,(x) ^ [5f J(x)]'^'^^W(x) , K^{x) ^ [Sy.ix)^ hlf\x) , 

(21) 



with 



[Sn{x)] = P exp (^tg dtn- Al{x + tn) tl)j , 



[Sf^{x)Y'' = Vexplig [ dtn- v4^(x + tn) {-1^^) 



oo 



[Sy, (x)] ''' = V exp (^-tg dt vs ■ A1{x + tvs) t^,^ , (22) 

and similarly for the ra-coUinear and h^^ fields. We have used the superscript "adj" to indicate 
Wilson lines in the adjoint representation. The fields with the superscript (0) no longer interact 
with soft gluon fields. Here are Gell-Mann matrices and /™'' are structure constants of QCD. 

One must supplement the effective Lagrangian with a set of operators describing the 
(qq, gg) — ti scattering processes. These operators appear in the effective Hamiltonian in 
convolutions along light-like directions with perturbative Wilson coefficients arising from in- 
tegrating out hard virtual fluctuations. We write this effective Hamiltonian as 



n,^{x) = J2f dt^dt^e^^^^^'^^^^y^ [cZ{hM)OZ{xMM) + CZ{tiM)OZ{xMM 

I,m 



(23) 



where / labels color structures and m labels Dirac structures. The operators can be written 
as 



0|l(x, ti, t,) = J2 (cf) X7 + hn) C xl' {x + t,n) Kl (x) ^ h 



{a} 



Ofmi^, tu t,) = ^ (cf )^^^ ^(x + t^n) AZ^{x + t,n) Kl{x) YZ K\{x) , (24) 

{a} 

where F^, F^, and F^ are combinations of Dirac matrices and the external vectors n, n, vs, 
and f4 (note that there can be contractions of Lorentz indices between F^ and F'4). The cj'^ 
and Cj^ are tensors in color space, whose indices {a} = {01,02,03,04} can be in either the 
fundamental or adjoint representation. For each channel, they are chosen to be the independent 
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color singlet (e.g. gauge-invariant) structures needed to describe the scattering amplitude. We 
will choose the color structures to be in the singlet-octet bases 

When squaring the amplitude and summing over colors, one must evaluate products of the 
color structures with their indices contracted. In the absence of soft gluon emissions, these are 

of the form icV^ (c^?) (and similarly for the gluon fusion channel) and are equal 

\ J / 01020304 \ -J / 01020304 \ ■J / ^ 

to an A^-dependent factor multiplying 5ij. In this sense the color structures are orthogonal, 
but not orthonormal. 

Time-ordered products of the operators (24) with the SCET and HQET Lagrangians de- 
scribe the coUinear and soft contributions to the [qq, gg) — )■ ttX scattering amplitudes in 
QCD, where the final state X contains any number of soft gluons from real emissions. In 
the formulation used so far, the final state X is built up through insertions of the interaction 
Lagrangian (20) into the SCET operators. To account explicitly for soft gluon emission to all 
orders in the strong coupling constant, it is convenient to use the decoupling relations (21) 
and represent the soft gluon interactions by Wilson lines. Performing this decoupling and 
dropping the superscripts on the new fields, the operators factorize into products of collinear, 
heavy-quark, and soft-gluon operators in the form of Wilson loops. The resulting operators 
in the qq channel read 

Oj^{x,t,M = E (^^W [Otix)]"'' mx,t,Mf'' mx)]^""^'^'^ , (26) 

{a},{b} 

where 

K(x)] '^'^ = h% (x) C (x) , [O'^ix, t,,t,)] = {x + hn) T'^ x'^ {x + hn) , 

[Qs^^)] w,w ^ [5;(^)]''3- [^^^(^)]-^* [slix)]"'"' [SMr"' . (27) 

Those in the gluon-fusion channel are obtained by making the obvious replacements in Dirac 
structure and the collinear operators, and by changing the Wilson lines Sn and Sn to the 
corresponding ones in the adjoint representation. 

In the form shown above, the different sectors no longer interact with each other. After 
squaring the amplitude, this property leads to the factorized form (10) for the hard-scattering 
kernels. A complication is that the color indices on the Wilson lines representing the soft-gluon 
interactions act on the color structures (25) and can mix them into each other. Consider, for 
example, the calculation of an 0{as) correction due to soft gluon exchange between partons 
1 and 3 in the squared amplitude in the qq channel. This could either be from the product 
of diagrams involving the real emission of one soft gluon, or from a virtual diagram with a 
soft loop (which would be scaleless for on-shell quarks, but appears in the calculation of the 
anomalous-dimension matrix). After summing over colors one must evaluate contractions of 
the form 

/ {bi 026304}^! "sbs V^J y {01020304} ■ ^ 
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In general, this contraction is not proportional to Sjj, and must be worked out case by case. 
This mixing of the different color structures due to soft gluon exchange is responsible for the 
non-trivial matrix structure of the hard and soft functions. To organize this color algebra, it 
is convenient to use the color-space formalism [77, 78] . Before moving on to the calculation of 
the differential cross section in Section 3.3, we briefly pause to review this formalism. 

3.2 Color-space formalism 

Consider the on-shell scattering amplitudes for {qq, gg) — )■ tt, for a given color configuration 
of the external particles. We write this in the quark channel as 

Mf^y = {t'HPs) t"^(P4) I ^eff(O) I q'^^ipi) r{P2)) , (29) 

and also define the object A^j^} in the obvious way. In what follows, we will drop the su- 
perscript indicating the channel, and work with a single amplitude which can represent either 
Ji^m QY _A/f99. As in [77,78], we introduce an orthonormal basis of vectors {|ai, 02, 03, 04)}, 
where the indices {a} refer to the colors of the external particles. The amplitude can then be 
written as 

-M{a} = (01,02,03,04 1 A^) , (30) 

where the object is an abstract vector in color space. Since we only consider color-singlet 
amplitudes, we can decompose the QCD amplitude into the set of color structures (25). In 
the color-space formalism, this is done by writing 

i-^) = E-^^EHwiw) = -^^|c^) ' (31) 

/ {a} 

where the coefficients Aij are combinations of Dirac matrices, external vectors, spinors, and 
polarization vectors. Note that with this definition the basis vectors \ci) are orthogonal but 
not normalized, so to project out the A^/ one must use 

Mi = -^{ci\M) . (32) 
The square of the amplitude summed over colors is then given by the inner product of |A1): 

colors / / {a} 

Following [77,78], we introduce color generators Tj to describe the color algebra associated 
with the emission of a soft gluon from parton i = 1,2,3,4. These matrices act on the color 
indices of the i-th parton as 

Tn...,a,,...) = {T^^),^aA...,k,...) . (34) 

If the i-th parton is a final-state quark or an initial-state anti-quark we set {Tf)i,a = t^^, for a 
final-state anti-quark or an initial-state quark we have {T[)ha = —t^b, and for a gluon we use 
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{T^)ba = if""^^- We also use the notation Tj-T,- = T^T^, and denotes the quadratic Casimir 
operator in the representation of the i-th parton, with eigenvalues Cp for quarks and Ca for 
gluons. Since we consider color-singlet amplitudes, color conservation implies the relation 

(Ti+T2 + T3 + T4)|7W) = 0. (35) 

We will be particularly interested in products of color generators acting on the scattering 
amplitudes, which appear in the calculation of perturbative corrections to the differential cross 
section. We write such products as, for instance, 

{M I T2 ■ T4 I M) = A^:,,,,,,,(T2'^)..a.(T4=)6.a.-Ma.a.a3a4 • (36) 

Rather than evaluating such expressions for each amplitude A^, it is more convenient to work 
out how the products ■ Tj act on the basis vectors \cj): 

T,-T,\cj) = [T,-T,]^j\cj) , (37) 

where the matrix elements are given by 

1 

(c/ I C/) 



[^^ ■ T^iJ = TTTTX (c/ I T^ ■ Tj I cj) . (38) 



It is worth emphasizing that while the generators themselves act in the abstract color space, 
on the left-hand side of the above equation they are just labels to identify a matrix acting in 
the space of color-singlet structures. This matrix is thus a 2 x 2 matrix for the qq channel, and 
a 3x3 matrix for the gg channel. 

We now consider the SCET representation of the amplitude Ai{a}- For the qq channel, 
this is equal to 

^ / dtidt2Cim{ti,t2) (C/)|^-^ 

t'^iP3)i'^{p.) M'^'^[O^,]'^'^[of'^^'^^^(0,ti,t2) q''^{pi)q'^{p2)) , (39) 



where no summation over the set of indices {b} is performed. We have made clear that for 
the amplitude to be non-zero, the colors of the heavy-quark and collinear fields must coincide 
with those of the external partons. The reason is that after the decoupling of soft gluons, the 
heavy-quark fields are effectively free fields, while collinear exchanges take place only within 

itself and are diagonal in color space. We can therefore suppress the color indices on 
and as well as on the external states, and keep in mind that we shall sum over the color 
indices of O'^ when we square the amplitude. Color correlations such as (36) are mediated by 
the exchange of soft gluons, which are represented by the Wilson lines in the soft operator Os- 
To describe these exchanges in the color-space formalism we use the operator 

Osix)= [S^SlSlS,,]ix) , (40) 

where the Wilson lines Si are defined as in (22), with the color generators promoted to the 
abstract ones Tj according to the rules stated below in (34). The decoupling relations (21) 
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are of the same form for quarks and gluons when expressed in terms of the Si [79], so this 
operator is used for both the qq and gg channels. Its action on the basis vectors |c/) is then 
defined according to (36). We define the full SCET operator as Om = O'^O'^O^ . 
To evaluate the partonic matrix elements of the coUinear fields we use 

(O I {XnY^itn) I Pu Si) = Saai e"**"'^ Ua{Pi, Si) , 

(0 I {An±)litn) I Pi; ai, s,) = 5aa, e-''^-P e^iPi, Si) . (41) 

The heavy-quark fields are always taken to be on-shell, so their partonic matrix elements are 
equal to HQET spinors multiplied by Kronecker delta symbols in the spin and color indices. 
Upon taking the partonic matrix elements, the integrals over ti and t2 produce the Fourier- 
transformed Wilson coefficients 

CUM,mt,cose,fi) = J t/tirft2e-^*^"-^^-**^"-P^ (7,^(^1,^2), (42) 

where we have made the full dependence of the momentum space Wilson coefficients on the 
kinematic variables explicit. Defining a vector of Wilson coefficients as 

\C^) = Y^Ci^\cj) , (43) 

and introducing the symbol ((. . .)) for partonic matrix elements as in [55], i.e., 

((0„)) = (tip-s) tip,) I 0^(0, 0, 0) I q{p,) q{p2)) , (44) 
we can write the color-space representation of the SCET scattering amplitude as 

\M) = $]((0„,)) \Cm) . (45) 



3.3 The differential cross section 



We now return to the derivation of the factorization formula (10) by calculating the partonic 
cross sections. The hadronic cross section is then obtained by convoluting these results with 
the PDFs. Below we will discuss the qq case in detail; the gg channel can be analyzed in an 
analogous way. 

The differential cross section is given by the phase-space integral of the squared amplitude 
^'^ = Y H2t)^2E, {2tf2E. +P3 +P4 -Pi -P2) 

2 

5^(t(p3)t(P4)^s(p.)|O^(0)|g(pi)g(p2))|CJ , (46) 



44 



where Om{x) = Om(a;,0,0), and we have used translational invariance to write the result 
in terms of the Fourier-transformed coefficients Cm from (42). The cross section is implicitly 
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summed over the external colors and spins, and the factor = N for quarks and dji = N"^ — 1 
for gluons arises from averaging over the colors of the initial-state partons. The factor 1/4 
accounts for the averaging over the polarizations of the initial-state partons. Since the different 
types of effective-theory fields do not interact with each other, we can factorize the matrix 
element into soft, heavy-quark, and coUinear pieces. The partonic matrix elements of the 
heavy-quark and collinear pieces just give the usual products of spinors, which we combine 
into the tree-level matrix element 

((0.„))t..cc = {{Ot{0)O';MKoe = (tips) Kp,) I Ot{0)O:^{0) I g(P2)>,,,, . (47) 

Summed over spins, these give rise to the usual Dirac traces. The matrix element of the soft 
operator is taken using the vacuum as the initial state and is of the form | {Xs{ps) \ 0^(0) | 0) p. 
This can be evaluated directly, but we prefer instead to sum over the final states Xs and convert 
it into a forward matrix element using the formalism explained in Appendix C of [64] . Written 
in this form, the differential cross section reads 

da = [ rf4^e^{Pi+P2-P3-P4).x 

2s (271)32^3 (27r)32E4 



4rf 



((0™))Lc ((an'))u-ec X (O | T[0^t(^)] T[O%0)] \ 0> 



(48) 



where T and T represent time and anti-time ordering [64]. Since we treat the soft scale as 
perturbative, the soft fields can be integrated out by evaluating the vacuum matrix element. 
We now define a hard matrix and a position-space soft function as 



H {M,mt, cos 9, fl) = | 7^ V ((Ow))trcc \Cm') {Cm\ {{Om))i 



tree 5 



m,m' 



Wix, /i) = ;^ (0 I T[0'H^)] T[O^(0)] I 0) , (49) 



R 



where we have chosen the prefactors to match the overall normalization of (10). The elements 
of these matrices in the chosen color basis are defined as 

Hij = , I I , {ci\H\ cj) , Wij = {ci\W\ cj) , (50) 

{Cl\Cl) {Cj\Cj) 

so that the term in square brackets in the second line of (48) is proportional to 

Tr[ifl¥] =Y,HijWji. (51) 
i,j 

In order to compute the invariant mass spectrum for the ti pair, we define q = p^ + Pi and 
insert 

1 = j d^q dM (g - P3 - Pi) 2M 8{M^ - q^) (52) 
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into (48). After performing the integral using the first (5-function, the integration usin^ 
the second, and carrying out the trivial angular integration in the ps integral, we arrive at 

(fa 16M 1 f d^q 



dE,\ps\ 6{M'-2q.p,) 



dMd cose 3s (27r)' 

X j d^x e^^fi+P^--^)-^ Tr [H{M, mt, cos 9, /i/) Wix, /i/)] , (53) 



where {p^l = — mf and go = ^/ M'^ + (p. In the partonic center-of-mass frame, we have 
|g| = (9(a/S(1 — z)), so we can set qo = M and drop q in the 5-function. Then the integral 
over q produces a factor 6^{x) from the exponential, and after a few manipulations we find 

Tr [ii'(M,mt, cos e,/i/) 5(71(1 - z),M,mt, cos e,/i/)] , (54) 



dMd cos 9 3sM 
where the momentum-space soft function is defined by [64] 



S {y/§{l -z),M,mt, cos 9, ^i) = ^/I [ ^e'v^^^"")""/' Vr(xo,f = 0,/i) 

J 47r 



(55) 



The hadronic cross section (5) is now obtained by convoluting the partonic cross sections 
in (53) with the parton luminosities. Comparing with (10), we finally arrive at the factorized 
form of the hard-scattering kernel 

C{z, M, mt, cos 9, fif) = Tr [H{M, mt, cos 9, fij) S{\/1{1 — z),M,mt, cos6', /i/)] , (56) 

where we have set s = everywhere except in the first argument of the soft function. 

4 The hard and soft functions at NLO 

In this section we describe the calculation of the soft and hard matrices up to NLO in pertur- 
bation theory. 

4.1 Hard functions 

The hard functions are related to products of Wilson coefficients, as shown in (49). To obtain 
the Wilson coefficients Cjm, one matches renormalized Green's functions in QCD with those 
in SCET. The matching can be done with any choice of external states and infrared (IR) 
regulators. It is by far simplest to use on-shell partonic states for {qq, gg) — )■ ti scattering and 
dimensional regularization in = 4 — 2e dimensions to regularize both the ultraviolet (UV) 
and IR divergences. With this choice, the loop graphs in SCET are scaleless and vanish, so 
the effective-theory matrix elements are equal to their tree-level expressions multiplied by a 
UV renormalization matrix Z. The matrix elements in QCD, on the other hand, are just the 
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virtual corrections to the (gg, gg) — )■ tt scattering amplitudes. The matching condition then 
reads [55, 74, 79] 

limZ-i(e,M,mt,cos^,/i) |7W(e, M, m^, cos ^)) = V((0„))trec |C^(M, m^, cos ^, /i)) , (57) 

e-5>0 ^ — ' 

m 

where A4. is the UV-renormalized virtual QCD amplitude expressed in terms of with rti = b 
active flavors. We have moved the SCET renormalization matrix Z to act on the QCD 
amplitude, so that both sides of the equation are finite in the limit e — )■ 0. The explicit results 
for the matrix elements Zu in our color basis for the qq and gg channels can be found in [55] . 

In practice, we are not interested in the Wilson coefficients themselves, but rather the hard 
matrix Hij. To calculate this, we first define 



|A^,en) = limZ-i(e) \M{e)) = 47ra, |A^W> + ^ jA^^) 
e-s>o L' 47r 



and expand the hard function as 



H = ai 



8df 



(58) 



(59) 



Using (57) and (58) to express the SCET matrix element in terms of the finite, IR-subtracted 
QCD amplitudes in the definition of the hard function (49), the matrix elements (50) can be 
written as 



H 



(0) 

ij 



1 



4 {ci\ci) {cj\cj) 



Cl 



Mi'i) {Mit I cj) 



H 



(1) 
IJ 



1 



{cj I Mi^l) {Mill I cj) + {cj I A^(^) > {Mi^l I cj) 



4 (c/lc/) {cj\cj) _ 

The leading-order result for the qq channel follows from a simple calculation and reads 



(60) 



H. 



(0) 
QQ 




M4 



M2 



(61) 



while that for the gg channel is 



/f(0) 

aa 



( ^ 



1 tj— Ml 

N AP 
J_ 

N 



1 fj— Hi 

N M2 



hp 



J 



2tiUi 



tl + ul 4mj 4m^ 



M4 



+ 



M2 tiMi 



(62) 



To calculate the NLO hard function requires the one-loop virtual corrections to the partonic 
scattering amplitudes, decomposed into the singlet-octet basis. Although results for the NLO 
virtual corrections interfered with the Born-level amplitudes exist in the literature [13-15], 
results for the one-loop amplitude decomposed into our color basis are not available and must 
be calculated from scratch. For this purpose we use in-house routines written in the computer 
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algebraic system FORM [80]. The results are rather lengthy (especially for the gg channel) 
and are in a Mathematica package which can be obtained from the authors upon request. 

We have been able to perform several checks on our results. First, we have verified that 
applying the renormalization factor Z to the tree-level amplitude indeed absorbs the IR poles 
in the UV-renormalized QCD amplitudes at one-loop order. Second, we have checked that 
inserting the results for the products of one-loop hard functions and tree-level soft functions, 
given in (65) below, into the formula for the differential cross section, we reproduce the results 
of [13-15], as required. Finally, using the one-loop hard functions, we were able to calculate 
the IR singularities of the two- loop amplitudes in [55] , which agree with all the available results 
in the literature [22-24]. 



4.2 Soft functions 

The soft functions are given by the vacuum expectation values of the soft Wilson-loop oper- 
ators, as defined in (49). In what follows we will calculate the one-loop corrections to these 
objects directly in position space. When performing the resummation in the next section, it 
will be more convenient to work with the Laplace-transformed functions. They are defined as 

s{L, MjiJit, cos9, fi) = --j= J du exp ^ — ^ ^ S{uj, M, nit, cos 9, fi) 



-2i 



(63) 



where the second equality was shown in [64] and follows from the functional form of position- 
space Wilson loops [81]. 

We expand the soft functions in power of as as 



5(0) 



Air' 



^(1) 



47r 



+ . . . . 



(64) 



At leading order, the Wilson loop is just the unit matrix, so sf] = (c/ | cj) /du, and it is easy 
to show that 



^(0) 






(N 





\ 


s(o) _ 





N 
2 







\o 





Ar2_4 
2N / 



(65) 



At NLO, the soft functions receive contributions from the diagrams depicted in Figure 1. The 
calculation is similar to that in [81]. To evaluate the diagrams we associate an eikonal factor 
v^/k-Vi multiplied by a color generator Tj for each attachment of a gluon to a particle with 
velocity Vi (we define vi = n and V2 = n), and contract with the cut gluon propagator in 
position space, which in Feynman gauge reads 



D 



X) = -g 



d'^k 
(27r)° 



-ik-x 



(271) 5{k^) eik"" 



(66) 



15 



^13 



Xl4 




^24 




Figure 1: Diagrams contributing to the soft functions at NLO. The thick hues represent Wilson 
hues in the time-hke directions and W4, the thin hues Wilson lines in the light-like directions 
n and n, and the cut curly lines represent the cut gluon propagator (66). 



We can then write the bare soft function in position space as 



(67) 



where the matrices Wij are related to products of color generators and will be given in (72) 
and (73) below. The integrals Xij are defined as 

-ik°xo 



Xij(e,Xo,/i) 



(47r/i 



2\e 



IT 



2-e '^J 



Vi ■ V,; I d k 



(27r) 5{e) e(fcO) 



(68) 



Vi - kvj ■ k 

which are obviously symmetric in the indices i and j. The integrals Xn = X22 = 0, and the 
non- vanishing integrals are 

I12 



vr' 
"6 



— H — Lq + Lq + — , , 



2 2 

T33 = X44 = — h 2Lq — — In , 
e Pt 
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+ r/2 \ 27?^ 

— ^ i- + 2Lo\ Inx, -In^x, + 41na;,ln(l - x,) +4Li2(xs) - — 

- Lo - n ^- + — + 2L 2 ^ + 2L 2 — — - 



2^14 — 2^23 — Xi^{yt — Zu) , 

where = (1 - + A), Z/t = -ti/rri^ - 1, Zu = -Ui/rn^ - 1, and 



Lq = In 



(69) 



(70) 



The renormahzed soft functions W^^^ can then be obtained by subtracting the divergent part 
from Wjl^^p. Later on we will need the Laplace-transformed function s, which according to 
(63) is obtained by replacing Lq —L. To finish the calculation, we must also determine the 
matrix elements of 



For the qq channel, the results are 



u = -T- (c/ I T^ ■ Tj I Cj) . 



"^12 — "^34 



Cf Un^ 

"4iV I -1 



(71) 



"^33 ~ "^44 



Cf 2N 
^ \ CfT 



QQ QQ 
'^13 = ^24 



Cf 



1 



2 \1 2Cf- -^ ' ' 



QQ QQ 

wll = w^l 



while for the gg channel we obtain 



C, 



2N 



w 



gg 

12 



w 



gg 

34 



1 


/4iV2 












4 






V 





/CfN 










1 

4 


\ 
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-iV' 
-A^ 1 



\ 



2_4 

\ 

/ 





Ar2_4 



4Af2 



(72) 



w 



99 
33 



W 



99 
44 



2N 



w 



99 
13 



W 



99 
24 



W 



99 
14 



W 



99 
23 



/2iV2 


V 

/ 

AN 

Vo 
/ 

-AN 
\ 





N^ 
_ 



4iV \ 

N^-A N^-Aj 



-AN 
-(iV2 _ 4) 



\ 

N^-A ) 



(73) 



The NLO contributions to the cross section from real emissions in the soft hmit have been 
known for some time [18]. Transforming our results to momentum space, we have checked 
that the integrals (68) are consistent with those given in Appendix A of that paper (after 
taking into account some misprints in [18] later corrected in [82]). As another check, we have 
used our results along with the one-loop hard functions from the previous section to calculate 
the total partonic cross sections at NLO in the limit /3 — > (see Section 6.3), reproducing the 
analytic expressions from [16]. Finally, using the RG invariance of the cross section we will 
derive the RG equation for s in the next section. We have checked that our one-loop result 
satisfies this equation, which also justifies our procedure of simply subtracting the 1/e poles 
in the bare function to get the renormalized results. 



5 Threshold resummation in SCET 

In the region where the cross section is dominated by the threshold terms, one needs to resum 
the leading singular terms in (1 — ^) to all orders in perturbation theory. This is accomplished 
by deriving and solving RG equations for the hard and soft functions in the effective theory, 
which will be described in what follows. Since these equations contain information on the 
logarithmic structure of the hard-scattering kernels at higher-orders in perturbation theory, 
they can also be used to derive an approximate NNLO formula for the differential cross section 
in the threshold region. We discuss this further in Section 5.2. 

5.1 RG evolution and resummation at NNLL 

The hard function satisfies the evolution equation 

— H(M, rrit, cos 9, fx) = Th(M, nit, cos 6, /i) H(M, nit, cos 9, fx) 
(iln/i 

+ H{M, mt, cos 9, /i) (M, mt, cos 9, /i) . (74) 

Using (50), we can write the above equation in a matrix form, where the matrix elements of 
Fff are defined according to (38), and can be obtained from the matrices Tgg or Tgg in [55]. 
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The form of the evolution equation follows from (58) and (60), along with the defining relation 



d 



H 



(75) 



for the anomalous dimension. The explicit results to two-loop order are 



qq 



Cf7cusp(«s) ( 



^\ ( ^^^ 



iTx + Ci.7cusp(/334, as) + 275(a,) + 27«(a,) 



0' 
1 



+ «7r -7cusp(/334,aJ 



f2 

+ 7cusp(as) In -\ 



^ 

1 



2 

-N 



^ 

" 2 



(76) 



and 
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iV7cusp(as) ( In 



M2 
7^ 



~2 



7cusp(as) In 



in]+CF 7cusp(/334, as) + 27^(a,) + 2^^ {as) 



/o o\ 
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vo i; 



+ «7r -7cusp(/334,as) 



+ 7cusp(as) In^ 
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(77) 



where the various anomalous- dimension functions can be found in the Appendix, and the cusp 
angle ^^.a = in — ln(l + (5t) /(I — Pt)- The solution to the evolution equation can be written as 

H{M, rrit, cos 6, /i) = U{M, mt, cos 6, /x/i, /i) H{M, nit, cos 6*, Hh) U\M, mt, cos 9, Hh, n) , (78) 

where the unitary matrix U satisfies the equation 



d 



(iln/i 



U{M, nit, cos 9, fih, /i) = Th{M, nit, cos 9, /i) U{M, nit, cos 9, fih, ^J) 



(79) 



The matching scale niust be chosen of order a typical hard scale, so that the matching 
condition for the hard function is free of large logarithms. With the help of the evolution 
matrix U , the hard function can then be evolved to an arbitrary scale yU.. The formal solution 
to this equation is 



/djj! 
— Th{M, rrit, cos 9, /i') 
/i' 



(80) 
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where the path-ordering is necessary because is a matrix. To evaluate the path-ordered 
exponential, it is convenient to separate the explicit logarithmic dependence on the scale n, 
which is related to Sudakov double logarithms, from the remaining piece, which is related to 
single logarithmic evolution. We thus write the anomalous dimension as 



TniM, mt, cos 9, /i) = Tcuspla^*) (in - in J 1 + 'y^{M, rrit, cos 9, as) 



^1) 



where Fcusp is equal to Ci?7cusp for qq and A^7cusp for 99, and the matrices 'y^ are defined 
through a comparison with (76) and (77). Since the term proportional to Fcusp multiplies the 
unit matrix, we can factor this piece out of the path-ordering and evaluate it using standard 
techniques. The result for the evolution matrix is then 



C/(M, nit, cos 9, fifi, jj) = exp 



2S'(/i/j, yu) - ar(/ih, fi) [\n — -iir 



M 



2 



u{M, rrit, cos 9, jj^h, jJ^) ■ 
(82) 

The RG exponents in the square brackets of the exponential factor are given by 



S{f^„f^) = - j d^^-'ff^ I ar(/..,/.) = - J da^^^ , (83) 

where P{as) = dag/dlnfi is the QCD /3-function. The quantity u contains the non-trivial 
matrix evolution due to 7'^ and reads 

/diCy. 
7^(M, mf, cos6', a) . (84) 
P(a) 

The perturbative solutions to the above equations are reviewed in the Appendix. 

We now turn to the evolution of the soft function. We derive its evolution equation by 
using the RG invariance of the cross section, 

TT[HS]0ff = O, (85) 



along with the evolution equations for the hard function and parton luminosities. The evo- 
lution equation for the hard function was given above, and the parton luminosity functions 
satisfy the DGLAP equations [83-85]. While the full DGLAP equations involve fiavor mixing, 
what we need here is the x — 1 limit of them, which is fiavor- diagonal and can be written as 

d dr 

#(y,/i) = 2 / _P(x)#(y/x,/i), (86) 



c? In /i Jy ^ 

where P{x) is given by 



P(x) = HJ=eM + 2j*{a.) S(l - x) ^ 
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The evolution for the momentum-space soft function is then 
d 



(iln/i 



S{uj, M, rrit, cos 6, /i) 



2rcusp(«s)ln— + nit, cosd, as 

/i 



S{u, M, nit, cos 9, fi) 



2rcusp(«s) In - + Y{M, mt,cos6,as 



S{uj, M, nit, cos 6, fj) 

I S{u', M, nit, cos 9, /i) — S{u, M, mt, cos 9, /i) 



^Tcusp i^^s 



dco' 



u — u' 
where we have defined 

Y{M, mt, cos 9, as) = 1^{M, nit, cos 9, a,) + 27^(a,) 1 . 



(88) 



(89) 



As in [64], the non-local evolution equation for the soft function can be turned into a local 
one by the Laplace transformation (63). The evolution equation for the Laplace-transformed 
function reads 

s m — — , M, mt, cos 9, jj, 



rfln/i 



M2 , 

rcusp(as) In — + 7 (^, "^t, cos 9, 



s [ In — —, M,mt,cos9, fi 



M2 

— s ( In — —, M, mt, cos9, /i 



M2 



rcusp(as)ln^ +Y{M, mt,cos9,as 



(90) 



This can be solved using the same methods as for the hard function. Transforming the results 
back to momentum space, we find 

S{u, M, mt, cos 9, fif) = VI exp [— 4S'(yUs, ///) + 4a^<#>(/is, ///)] 
X u^(M, mt, cos9,jjLf, fis) s{dr^, M, mt, cos 6', /i^) u{M, mt, cos 9, jij, fig 



1 ( uj\'^'^ e-^T'^'' 



uj \iisj r(2?7) 

(91) 

where rj = 2ar(/is, /^/)- The soft scale fis should be chosen such that the contribution from 
the soft function to the cross section is perturbatively well-behaved, and will be discussed in 
detail in Section 6. 

Combining the results for the hard and soft functions, our final resummed expression for 
the hard-scattering kernel is 

C{z, M,mt,cos9, fif) = exp [4a^0(/is, /i/)] 



X Tr 



U{M, mt, cos 9, fih, fi.s) H{M, mt, cos 9, fih) U\M, mt, cos 9, fih, fis) 
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For values Hs < the parameter 77 < 0, and one must use a subtraction at 2 = 1 and 
analytic continuation to express integrals over z in terms of star (or plus) distributions [86]. 
Formula (92) can be evaluated order-by-order in RG-improved perturbation theory, using the 
standard counting In fih/f^s ~ ln(l — z) ~ l/a^. The perturbative solutions for the RG factors 
needed to evaluate the evolution matrix U to NLO in this counting scheme are given in 
(A. 2), (A. 3), and (A. 5) of the Appendix. The correspondence between this counting and the 
standard counting of logarithms (e.g. NLL, NNLL), along with the accuracy of the anomalous 
dimensions and matching functions needed at a given order, can be summarized as follows: 



RG-improved PT 


log accuracy 


r 

cusp 




H, s 


LO 


NLL 


2-loop 


1-loop 


tree-level 


NLO 


NNLL 


3-loop 


2-loop 


1-loop 



In the remainder of the paper we will use the logarithmic counting (e.g. NNLL) when referring 
to the resummed results obtained in this section. These results are valid for the leading- 
order term in the threshold expansion in (1 — z), whereas the full result at NLO in fixed- 
order perturbation theory also contains information on subleading terms. In phenomenological 
applications we can match the resummed results with the NLO fixed-order results to achieve 
an NLO+NNLL precision. The method for doing this is described in Section 6. 



5.2 Approximate NNLO results 



In the previous subsection we derived a formula for the resummed differential cross section, 
which is valid up to NNLL order. Starting from (92), it is also possible to obtain expressions 
for the differential cross section which are valid in fixed-order perturbation theory [56] . Indeed, 
our results allow one to obtain analytic expression for all of the coefficients multiplying singular 
plus distributions in the variable {1 — z) appearing in the hard-scattering kernels up to NNLO. 
With the same method, which is outlined below, it is also possible to determine analytically, 
up to 0{ag), the scale-dependent parts of the coefficient multiplying 5(1 — z). 

In order to derive fixed-order formulas from (92), we first set yU/i = /is = /i/ = /i. In that 
case the evolution matrix U is equal to unity, and 77 = 2ar(/i/,/is) — 0. The formula for the 
hard-scattering kernels then becomes 



C{z, M, rrit, cos 6, fi) = c{dri, M, rrit, cos 6, fi) i — ) 



12 J T{27]) (1 - z)l-2'? 



r)=0 



(93) 



where 



c(L, M, nit, cos 9, fi) = Tr H{M, nit, cos 9, fi) s{L, M, nit, cos 9, fi) 



(94) 
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By using (74) and (90) in combination with the analytic expressions for the hard and soft 
functions at NLO, it is possible to determine all terms proportional to In/i in the two-loop 
hard function H^'^\M, rrit^ cos 6', /i), as well as all terms proportional to L in the two-loop soft 
function s^'^\L,M, rrif, cos 0,/i). This information allows us to derive an approximate expres- 
sion for c at NNLO. By inserting that formula for c into (93), we obtain the corresponding 
NNLO expression for the hard-scattering kernel C. The results are conventionally written in 
terms of the plus distributions 



PJz) 



ln"(l - z) 
1-z 



(95) 



However, the right-hand side of (93) is more conveniently expressed in terms of the distribu- 
tions 



1 .^^nfM\i-zr 



1-Z 



IJ?Z 



(96) 



It is possible to show that taking the derivatives with respect to rj and the limit — in (93) 
is equivalent to making the following set of replacements in c(L, M, mt, cos^, /i): 

1 ^5(1-2), 

L^2Pl,{z) + 6{l-z) In — 

\ ^^ 



L"" ^ 4:P[{z) + 5{l - z) In^ 

GP^iz) - AiT^Poiz) + 6{l-z 
8P3(^) - 167r^P[{z) + 128(3^0(2) + 5(1 - z) 



In' 



^)+4C3 



In^ 



M2 

7^ 



+ 16C3ln(^ 



(97) 



In order to translate the P^ into the conventional P„ distributions, we employ the general 
relation 



fc=0 

fc-1 



1 — X 



The final result for the hard-scattering kernels at NNLO can be written as 
C^'^\z,M,mt, cose, fx) = D3 



'ln^(l - z)' 


+ D2 


'\n\l - z)' 






1-z 


+ 


1-z 



(9^ 
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ln(l - z) 



Dn 



+ CoS{l-z)+R{z) 



(99) 



The coefficients Dq, . . . , D3 and Cq are functions of the variables M, rut, cos 6*, and /i. The an- 
alytic expression for Di {i = 1,2,3) were first derived in [52] starting from resummed formulas 
in Mellin moment space. In [56] the coefficients Dq, . . . , D3 were completely determined by 
following the procedure outlined above. Dq can be calculated in this way because the process- 
dependent anomalous-dimension matrices in (76) and (77) are now known up to NNLO. With 
the same method, it was possible to calculate the scale dependence of (5-function coefficient 
Cq. The function R{z) is finite for z — >■ 1. The computation of the scale- independent part 
of Co requires the knowledge of the hard and soft functions at two-loop order. As long as 
these are missing, there is an ambiguity in Co, associated with the fact that it is possible to 
normalize the scale-dependent logarithms in an arbitrary way, i.e. 
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(100) 



Scale-independent terms proportional to the second logarithm on the right-hand side can be 
absorbed into the unknown //-independent part of the coefficient function Cq. Therefore, our 
numerical results for Cq depend on the choice of the second mass scale, which appears in the 
scale-dependent logarithms; we indicate this second scale by fiQ. In the next section we shall 
consider two different choices for /iq: fiQ = M (scheme A), and fiQ = rrit (scheme B). The 
situation is summarized in the following formula: 



r 4 



M 



In" — (scheme A) 



Co= I 



1=0 

4 



J2 Cn In" ^ (scheme B) 



i=0 



(101) 



The coefficients and (n = 1, 2, 3) are known, while and Cg are unknown. In [56], the 
explicit expressions for the coefficients Do, . . . and Co were collected in a Mathematica 
ffie, which can be downloaded from the arXiv version of that work. 



6 Phenomenological applications 

In this section we perform a numerical study of the invariant mass distribution, forward- 
backward asymmetry, and total cross section at the Tevatron and LHC The main purpose is 
to investigate the impact of the NNLL resummation and the approximate NNLO corrections 
on the central values and perturbative uncertainties in these observables. The higher-order 
corrections computed in this paper are limited to leading order in the threshold expansion, 
whereas the exact NLO results in fixed order also contain subleading terms in {1 — z). To make 
optimal use of our results, we match them onto the NLO fixed-order expressions in such a 
way that these subleading corrections are fully taken into account. For the resummed results. 
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NLO+NNLL accuracy can be achieved by evaluating differential cross sections according to 



^^NLO+NNLL ^ ^^NNLL 



I ^^NLO, sublcading 



= da 



NNLL 



da 



NLO 



leading 



(102) 



where da^^'^ is the exact result in fixed order, and (icr^'"°'^°^'^™s|^^ = da^^^^\^^=^j_^=^j_j captures 
the leading singular terms in the threshold limit z — > 1 at NLO. The term da^^^' subieadmg 
of subleading order in {1 — z) and ensures that the total result reduces to the exact fixed-order 
result when all the scales are set equal. It also makes the result invariant under variations of 
the factorization scale, up to terms at NNLO in the perturbative expansion, even at subleading 
order in (1 — 2;). To obtain approximate NNLO results in fixed order, we simply add the NNLO 
correction onto the exact NLO results, i.e., 

^^NNLO, approx ^ ^^NLO ^ ^^(2), approx ^ (^^gg^, 

where da^"^^' '^pp™^ is the NNLO correction to the differential cross section obtained using the 
coefficient function (99), in either scheme A or B described at the end of the previous section. 

It is important to examine the relative size of the leading and subleading terms in (102). If 
the subleading terms were comparable in size to the leading ones, it would make little sense to 
resum the logarithms in (1 — z), which we have discussed so far, or to construct approximate 
fixed-order formulas at NNLO which capture only the effects of the singular terms. The naive 
expectation would be that the singular terms are dominant only when r = M^/s — )■ 1, since 
then the integrand in (5) is needed only in the 2; — )■ 1 limit and the less singular terms are clearly 
subleading. However, the most interesting region for phenomenology ranges from M ~ 2mt 
to around 1 TeV at the Tevatron and up to several TeV at the LHC, which corresponds to 
r < 0.3 (at most). For the leading-order singular terms to be dominant, it is necessary that 
the parton luminosity functions ffij^r/z,^) fall off sufficiently fast for t/z — )■ 1 that only 
the largest values of z give significant contributions to the integrand, an effect referred to 
in [64-66] as dynamical threshold enhancement. The study of the invariant mass spectrum 
in [56] indicated that such an enhancement actually does take place in the range of invariant 
mass from near the peak of the distribution at M ~ 380 GeV up to much higher values. 
This is illustrated in Figure 2, which shows the invariant mass distributions at the Tevatron 
and LHC predicted using different approximations in fixed-order perturbation theory. The 
difference between the boundaries of the dark NLO bands and the dashed lines is due to the 
small contributions from the subleading terms c?crNLO,subicading (102). The fact that, even 
at these relatively low values of M, the leading terms provide a very good approximation 
to the full NLO result provides a strong motivation to study within our formalism higher- 
order corrections to integrated quantities such as the total cross section and forward-backward 
asymmetry, which receive their dominant contributions from low values of the invariant mass. 

We will always do the matching onto fixed-order results as in (102) and (103), when the goal 
is to provide quantitative phenomenological predictions. Such a matching is straightforward 
for integrated quantities such as the total cross section and forward-backward asymmetry, 
since the NLO results in fixed order are available in analytic form. For the invariant mass 
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Figure 2: Fixed-order predictions for the invariant mass spectrum at LO (light bands) and 
NLO (dark bands) for the Tevatron (left) and LHC (right). We use MSTW2008NLO PDFs [87] 
with as{Mz) = 0.120. The width of the bands reflects the uncertainty of the spectrum under 
variations of the matching and factorization scales. The dashed lines refer to the leading terms 
in the threshold expansion. 



distribution, on the other hand, the fixed-order NLO results are available in the form of Monte 
Carlo programs such as MCFM [88]. This makes it difficult to get accurate values of the top- 
quark pair invariant mass spectrum at high M, where the differential cross section is small, 
and makes it impractical to calculate the spectrum with the scale choice fif = M used in 
the next section, since doing so would require to run the program separately at each point 
in fif. (Monte Carlo programs generate the invariant mass spectrum by first producing a set 
of events for a given fif, and then grouping them into bins in M). When we study certain 
aspects of the invariant mass distribution in Section 6.1, we will take the NLO correction in the 
threshold approximation, so that (102) and (103) are evaluated with da^^'~' — )■ da^^^' ioadmg_ 
This is still a good approximation to the full NLO result, and allows us to study the qualitative 
behavior of the invariant mass spectrum with Hf = M over a large range of M, as well as 
PDF uncertainties, in a simple way. For this purpose, we also define an NNLO approximation 
which includes only the singular terms at threshold in the NLO correction: 

^^NNLO, leading ^ ^^NLO, leading ^ ^^(2), approx _ /^q^n 



6.1 Invariant mass distribution: Systematic studies 

The invariant mass distribution is obtained by integrating the doubly differential rate over the 
range —1 < cos6' < 1. The resummed results (92) depend on the three scales ^s,fJ'h, and fj,f, 
and to give a numerical result we must first specify how to choose them. In the similar cases 
of Drell-Yan [64] and Higgs production [65, 66] at threshold, the soft and hard scales were 
chosen by examining the contributions of the one-loop soft and hard matching coefficients as 
functions of the scales fis and fih, and then choosing default values of the scales in such a way 
as to minimize these corrections. We shall use this approach here, a small complication being 
the extra dependence on the kinematic variable cos^ in (92). For the analysis of this section 
we use the default set of MSTW2008NNLO PDFs [87], take a,(Mz) = 0.117 with three-loop 
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Figure 3: Top left: Ratio of the one-loop correction from the soft function over the leading- 
order result for top pair production at the Tevatron, as a function of /i^/M, for M = 400 GeV 
(dark), M = 700 GeV (medium), and M = 1000 GeV (light). Top right: The scale /i,/M 
determined by the point where the one-loop correction from the soft function is minimal, 
as a function of the invariant mass M. Bottom: Analogous plots for the LHC, but with 
M = 400 GeV (dark), M = 1000 GeV (medium), and M = 2000 GeV (light). 

running in the MS scheme with five active flavors, and employ the value rrit = 173.1 GeV for 
the top-quark mass defined in the pole scheme. Using a fixed set of PDFs helps to elucidate 
more clearly the behavior of the perturbative expansion of the hard-scattering kernels in higher 
orders of perturbation theory. When presenting our phenomenological results in Sections 6.2- 
6.5, we will change the PDF sets according to the order of perturbation theory employed, and 
we will also study the theoretical uncertainties related to the parameterization of the PDFs. 
Finally, in Section 6.6, we will explore different schemes for the definition of the top-quark 
mass and study their impact on the phenomenological results. 

6.1.1 Determination of the matching and factorization scales 

We begin by examining the corrections from the NLO soft matching coefficient as a function of 
fig. We isolate this contribution by picking out the piece of the NNLL approximation to (92) 
proportional to s^^\ evaluating the differential cross section using only this piece, and dividing 
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Figure 4: Left: Ratio of the one-loop correction the from hard function over the leading- 
order result at the Tevatron, as a function of /i^/M, at M = 400 GeV (dark), M = 700 GeV 
(medium), and M = 1000 GeV (hght). Right: Analogous plot for the LHC, with M = 400 GeV 
(dark), M = 1000 GeV (medium), and M = 2000 GeV (light). 



the result by that at NLL, for the choice fij = fih = M. The results are shown in the left-hand 
plots of Figure 3 for the Tevatron and the LHC with -^i = 7 TeV, for several different values 
of M. We note that the corrections are larger at the LHC than at the Tevatron, especially at 
high values of M. This behavior appears to be a property of the gluon channel, which gives 
the dominant contribution at the LHC. The correction is generally at its minimum between 
M/4 and M/10, and moves to lower values of fig at higher values of M. The exact position 
of the minimum as a function of M is shown in the right-hand plots of Figure 3. To a good 
approximation, the numerical results for fig can be fitted by the function 

- (a + 6rV4)c l^U5) 

with a = —33, b = 150, and c = 0.46 for the Tevatron, and a = —1.3, b = 23, and c = 0.98 for 
the LHC at ^/s = 7 TeV. In Section 6.3 we will also study the total cross section at the LHC 
with y/s = 10, 14 TeV, and in those cases we use a = 0.95, b = 6.7, and c = 1.6. 

The most appropriate choice of the hard scale fih is not immediately apparent, since the 
invariant mass spectrum depends on the two hard scales rrit and M. As a guide to an appro- 
priate choice we look at the size of the correction from the hard matching function for different 
choices of /i^. We show in Figure 4 the correction obtained by isolating the contribution of 
H^^'' to the differential cross section at NNLL, and dividing it by the NLL result, for the 
choice fif = M and fig determined according to (105). We see that at lower values of fih closer 
to rrit the correction typically gets smaller and can even become negative. In this lower range 
of /i/i, however, the correction depends very strongly on the scale. The results are more stable 
in the range M/2 < fi^ < 2M, where the correction is generally below 30% at the Tevatron 
and between 20 — 40% at the LHC. In what follows we shall choose fi^ = M by default, in 
order to avoid the instability at lower fih. In the case of Higgs production, a negative hard 
scale squared /i^ ~ ''I^'h ~ chosen to minimize the logarithms arising from time-like 

kinematics [65,66]. In the tt case, however, there are both time-like and space-like momentum 
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Figure 5: Dependence of da/dM at the Tevatron (left) and LHC (right) on the scale Hf 
in fixed-order perturbation theory. The dashed-dotted lines show da^'^, the dashed lines 
^^NLO, leading^ ^^iq dark (light) solid lines the approximate threshold expansion (104) at 
NNLO in scheme A (scheme B). 



invariants, and it is not straightforward to tell which point in the complex plane should be 
chosen to minimize the logarithms. We have thus investigated the choice ^\ = M'^e^'^'^ with 
(ph varied between — vr and it. The results show that at the LHC, the correction is smallest 
for (ph ~ — TT, and is about 10% compared to 20% at (ph ~ 0. At the Tevatron, however, the 
minimal correction is obtained for values close to (ph ~ 0. In view of this, and since the cor- 
rections are in any case not very large, we will not go into this complication in our numerical 
analyses below. 

Finally, we must choose a default value for the factorization scale in both the resummed 
and fixed-order results. To do so, we study the behavior of the cross sections as a function 
of this scale. For the fixed-order results, the invariant mass spectrum as a function of /ij at 
M = 400 GeV and M = 1 TeV is shown in Figure 5. For the moment we do not match the 
results onto fixed-order perturbation theory at NLO, using instead the threshold expansion 
^^NLO, leading ^^NNLO, leading^ j^^ ^ result, our predictions are not strictly independent of the 
scale fif, but a slight scale dependence enters via subleading terms in (1 — z). At M = 400 GeV 
the approximate NNLO formulas differ from each other less at /i/- ~ 400 GeV than at ~ m^. 
The same is true at M = 1 TeV, but in this case the results become very unstable at yU/ ~ m^. 



29 



44 

^ 42 
O 40 

38 

^ 34 
b 

-XS 32 





M = 400 GeV 



















30 

100 200 



300 400 500 600 

tif [GeV] 



700 800 




500 1000 1500 



2000 




400 500 600 

/i/. [GcV] 



60 90 120 

^is [GeV] 



150 



Figure 6: Dependence of da/dM at the Tevatron on the scales fih, f^s, and fif in resummed 
perturbation theory. The defauh choices are Hh = IJ^f = M, and /is according to (105). The 
dashed-dotted hues refer to NLL, the dashed to NNLL. 



It therefore seems more appropriate to make the choice Hf ^ M when studying the invariant 
mass spectrum. The resummed resuhs at M = 400 GeV and M = ITeV as a function of /x/, 
with iih = M and fig as in (105), are shown in the upper two plots of Figure 6, for the case of 
the Tevatron (plots for the LHC would look very similar). Again, the results at /ij ~ M are 
more stable than at /ij ~ rrit, although compared to the fixed-order results the difference is 
less pronounced. We will thus make the choice /i/ = M by default in the resummed result. 

Having chosen default values for the scales, we now discuss in more detail the behavior 
of the fixed-order and resummed predictions for the invariant mass distribution. We have 
already seen how the results depend on the scale /ij. In the lower two plots of Figure 6, we 
show the dependence of the resummed results at M = 400 GeV on the scales fih and /i^. The 
results as a function of these two scales are significantly more stable at NNLL than at NLL. 

6.1.2 Convergence of the perturbation series 

An interesting difference between the fixed-order and resummed results is that the perturbative 
uncertainties and the size of the higher-order corrections in the resummed results depend much 
less on the value of M. This is seen in Figure 7, where we show the K factors and uncertainties 
in the invariant mass spectrum as a function of M, comparing fixed-order results with the 
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Figure 7: K factors {da / dM) / {da^'^''^'^^ / dM) in fixed-order perturbation theory (left) and 
{da/dM)/{da^^^''^^ydM) in resummed perturbation theory (right), for the Tevatron (top) 
and LHC (bottom). The hght bands in the fixed-order (resummed) results show LO (NLL) 
results, the medium bands show NLO leading (NNLL) results, and the dark bands in the 
fixed-order results refer to the approximate threshold expansion (104) at NNLO in scheme A. 
The width of the bands reflects the uncertainties associated with variations of the scales, as 
described in the text. 



resummed ones. The bands in fixed order reflect the uncertainty associated with varying 
the factorization scale around its default value /i/ = M by a factor of two. Here and in the 
following figures, for the approximate NNLO formulas we show results only in scheme A; those 
in scheme B look very similar. To make the bands in resummed perturbation theory, we first 
obtain uncertainties associated with the scales fif, fih, and fis by varying them individually up 
and down by a factor of two from their default values at each point in M, with the other two 
scales held fixed. We then obtain a total error by adding the three uncertainties obtained this 
way in quadrature. At the Tevatron, the K factors in resummed perturbation theory have 
smaller uncertainties and depend only weakly on M, compared to fixed order. The same is 
true at the LHC, although at small M the approximate NNLO results and the resummed ones 
have comparable uncertainties. 

To illustrate more precisely the quantitative differences between the various perturbative 
approximations to the invariant mass spectrum, we show in Table 1 the exact numerical values 
of the spectrum at the points M = 400 GeV and M = 1000 GeV. We have assigned uncertain- 
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Tevatron 


M = 400 GeV [fb/GeV] 


M= ITeV [fb/GeV] 
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[X [X Q+78+20 
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' ' -36-29 


7 /I q + 1.61 +0.69 
' -^^-1.07-0.70 

7 1 7+1.19+0.69 
-1.10-0.69 

1 n fiq+0.84 + 1.01 
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Q qQ+0.82+0.90 
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Q 1 1 +0.63+0.86 
-^•^-■--0.78-0.86 



Table 1: Values for da/dM for M = 400 GeV and M = 1 TeV at the Tevatron and the LHC 
at ^/s = 7 TeV. The first error refers to perturbative scale uncertainties, the second to PDF 
uncertainties, see text for a detailed explanation. 



Order 


PDF set 




LO 

NLO, NLL 

NNLO approx, NLO+NNLL 


MSTW2008LO 
MSTW2008NLO 
MSTW2008NNLO 


0.139 
0.120 
0.117 



Table 2: Order of the PDFs [87] and the corresponding values of the strong coupling used for 
the different perturbative approximations. 

ties associated with variations of the various scales by factors of two up and down from their 
default values. To obtain a total scale uncertainty for the resummed results, we have added 
the uncertainties associated with variations of fih, f^s, and /i/ in quadrature. We have also 
included uncertainties associated with the PDFs, by using the set of MSTW2008NNLO PDFs 
from [87] at 90% confidence level (CL). The perturbative scale uncertainties are smaller or 
comparable than those from the PDFs only once the NNLL or approximate NNLO corrections 
are taken into account. For the practical reasons explained earlier, we have not matched the 
higher-order results with the fixed-order NLO results. However, the threshold approximation 
works rather well. For reference, at the Tevatron the exact NLO results are {38.6~^l'l) fb/GeV 
for M = 400 GeV and (24.8^^8) " 10-3fb/GeV for M = 1000 GeV, while at the LHC they are 
(654+^^) fb/GeV for M = 400 GeV and (6.84+};^?) fb/GeV for M = 1000 GeV. The deviations 
from the leading NLO terms shown in the second line in both parts of the table are smaller 
than 7% for the Tevatron and 5% for the LHC. 
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Figure 8: Left: Fixed-order predictions for the invariant mass spectrum at LO (light), NLO 
(darker), and approximate NNLO (dark bands) for the Tevatron (top) and LHC (bottom). 
Right: Corresponding predictions at NLL (hght) and NLO+NNLL (darker bands) in re- 
summed perturbation theory. The width of the bands reflects the uncertainty of the spectrum 
under variations of the matching and factorization scales, as explained in the text. 



6.2 Invariant mass distribution: Phenomenological results 

After these systematic studies, we now present our final results for the ti invariant mass 
distributions at the Tevatron and LHC. Here and below, we will use different sets of PDFs, 
as appropriate for the order of the perturbative approximation employed. Strictly speaking, 
theoretical predictions obtained in resummed perturbation theory would require PDF sets 
extracted from data using resummed predictions for the relevant cross sections; however, such 
PDF sets do not exist at present. Since our resummed expressions include the bulk of the 
perturbative corrections appearing one order higher in Og, we use NLO parton densities for 
the NLO and NLL approximations, and NNLO parton densities for the approximate NNLO 
and matched NLO-I-NNLL approximations, as summarized in Table 2. The associated running 
couplings as{n) are taken in the MS scheme with five active flavors, using one- loop running 
at LO, two-loop running at NLO, and three-loop running at NNLO. 

We begin by studying in more detail the invariant mass spectrum at relatively low values 
of M, where it is the largest, in fixed-order and resummed perturbation theory. Contrary to 
the previous section, we now match the results in resummed perturbation theory with the 
exact fixed-order results at NLO using the MCFM program, according to (102). In this way 



33 



2.5r 



2.0 



Lo ^=1.96TeV 

NLO, leading 
NNLO, leading 



^ 1.5 




2.5 



2.0 



li:^ 1.5 



1.0 



400 500 600 700 800 900 1000 

M [GeV] 



0.5 



NLL 

NNLL 



^/s = 1.96 TeV 




400 500 600 700 800 900 1000 

M [GeV] 




> 

o 



100 



10 



b 

0.1 



0.01 



400 500 600 700 800 900 1000 

M [GeV] 



0.01 





v^= 1.96 TeV 


X 

NLL 




NNLL 


^x 
^x 
^x 
^x 
^x 
^x 
^x 
~~x 



400 500 600 700 800 900 1000 

M [GeV] 



Figure 9: Left: Fixed-order predictions for the K factor and invariant mass spectrum at 
LO (light), NLO (darker), and approximate NNLO (dark bands) for the Tevatron. Right: 
Corresponding predictions at NLL (hght) and NLO+NNLL (darker bands) in resummed per- 
turbation theory. The width of the bands reflects the uncertainty of the spectrum under 
variations of the matching and factorization scales, as explained in the text. 



we obtain state-of-the-art predictions, which include everything known about the perturbative 
series for the spectrum. Our results are shown in Figure 8. The bands reflect uncertainties in 
scale variations according to the same procedure explained in the previous paragraph, but in 
this case with /i/ = 400 GeV by default. For the range of M in the plot, this choice is very 
close to our preferred scheme /i/ = M, but allows for a simple matching with the fixed-order 
results from MCFM. One sees that the perturbative uncertainty estimated by scale variations 
is by far the smallest at NLO+NNLL order. 

We now consider the region of higher invariant masses, for which the dominance of the 
threshold terms is even more pronounced, as indicated by the convergence of the dark bands 
and dashed lines in Figure 2 toward higher M values. Figure 9 shows our results for the 
Tevatron, both in fixed-order and resummed perturbation theory. Figure 10 shows the cor- 
responding results for the LHC. It is impractical to match onto fixed-order results obtained 
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Figure 10: Left: Fixed-order predictions for the K factor and invariant mass spectrum at LO 
(light), NLO (darker), and approximate NNLO (dark bands) for the LHC. Right: Correspond- 
ing predictions at NLL (hght) and NLO-f-NNLL (darker bands) in resummed perturbation 
theory. The width of the bands reflects the uncertainty of the spectrum under variations of 
the matching and factorization scales, as explained in the text. 



using the MCFM program in this case; however, the differences compared with the shown 
curves are so small that they would hardly be visible on the scales of the plots. The upper 
two plots show K factors, which are defined as the ratio of the cross section to the default 
lowest-order prediction da^'^''^'^^ / dM . Contrary to Figure 7, we now use the same normaliza- 
tion in both fixed-order and resummed perturbation theory, so that the two spectra can more 
readily be compared to each other. The lower plots show the corresponding spectra directly. 
We observe similar behavior as in the low-mass region. The bands obtained in fixed-order 
perturbation theory become narrower in higher orders and overlap. The bands obtained in 
resummed perturbation theory are narrower than the corresponding ones at fixed order. The 
leading-order resummed prediction is already close to the final result. 

The information contained in Figures 8-10 can be represented differently in terms of the 
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Figure 11: Distributions da/dPt at the Tevatron (left) and LHC (right). 
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Figure 12: Comparison of the RG-improved predictions for the invariant mass spectrum with 
CDF data [9]. The value rrit = 173.1 GeV has been used. No fit to the data has been performed. 



very useful distribution da/d^t, with Pt defined as in (4). A simple change of variables yields 

do- _ 2mt(3t da 

dpr {I - dM ■ ^'''^ 

The resulting spectra for the Tevatron and LHC, obtained using RG-improved perturbation 
theory, are shown in Figure 11. As before, the distributions are normalized such that the area 
under the curves corresponds to the total cross section. Recall that the physical meaning of 
the variable /3t is that of the 3- velocity of the top quarks in the ti rest frame. The distributions 
show that the dominant contributions to the cross section arise from the region of relativistic 
top quarks, with velocities of order 0.4-0.8 at the Tevatron and 0.5-0.9 at the LHC. We will 
come back to the significance of this observation in the next section. 

In Figure 12, we compare our RG-improved prediction for the invariant mass spectrum 
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to a measurement of the CDF collaboration obtained using the "lepton + jets" decay mode 
of the top quark [9]. We observe an overall good agreement between our prediction and the 
measurement, especially for higher values of M. Apparently, there is no evidence of non- 
standard resonances in the spectrum. The only small deviation from our prediction concerns 
the peak region of the distribution, shown in more details in the right plot. This deviation 
has also been observed in [9], where a Monte Carlo study of the SM expectation has been 
performed. 

6.3 Total cross section: Phenomenological results 

The total cross section is obtained in our approach by integrating numerically the doubly 
differential cross section in the ranges — 1 < cos 6' < 1 and 2mt < M < y/s. In this case it 
is a simple matter to match onto NLO in fixed-order perturbation theory, using the analytic 
results of [16]. To do this, however, we can no longer correlate the factorization scale /ij 
with M, as we did when studying the invariant mass spectrum. Instead, we should resort to 
representative average values of M, which characterize the spectrum in the region yielding 
sizable contributions to the total cross section. One possibility is to take the location of the 
peak in the da/dM distributions, which is Mpeak ~ 375 GeV for the Tevatron and Mpeak ~ 
388 GeV for the LHC (see Figure 8). Another possibility is to take the average value (M) of 
the distributions, for which we find (M) ^ 445 GeV for the Tevatron and (M) ^ 496 GeV 
for the LHC. As previously, we take the fixed value /^/ = 400 GeV as our default choice. On 
the other hand, we are still free to choose the hard and soft scales as we have done so far 
and match with the fixed-order result as shown in (102). We display in Table 3 the central 
values and scale uncertainties for the total cross section obtained using this procedure. The 
results in resummed perturbation theory use Hh = M and /i^ chosen according to (105) by 
default, and the uncertainties are obtained by varying these scales and the factorization scale 
/// up and down by a factor of two and adding the different uncertainties in quadrature. The 
perturbative uncertainties in the fixed-order results are obtained by varying the factorization 
scale up and down by a factor of two from its default value. In addition to the perturbative 
uncertainties, we also list the PDF uncertainties obtained by evaluating the cross section with 
the appropriate set of MSTW2008 PDFs at 90% CL. As shown in Table 2, the LO cross 
sections are evaluated using LO PDF sets, the NLL and NLO cross sections using NLO PDF 
sets, and the NNLL and approximate NNLO cross sections using NNLO PDF sets. In the 
following tables, these different classes of predictions are separated by horizontal lines. 

A few comments are in order concerning the results shown in the table. At NLO the 
cross sections o"nlo. leading evaluated using only the leading singular terms from the threshold 
expansion reproduce between 95% (for the Tevatron) to almost 100% (for the LHC) of the 
exact fixed-order result at the default values of the factorization scale. The subleading terms 
in (1 — z), obtained by integrating cicTNLO, subicading, contribute the remaining few percent. In 
other words, the singular terms capture about 85% of the NLO correction at the Tevatron and 
practically 100% of it at the LHC. We cannot say whether the threshold expansion works so 
well also at higher orders in perturbation theory, although this does not seem unreasonable. 
Our best prediction is obtained by matching the fixed-order result with the resummed result 
at NLO-I-NNLL accuracy and is highlight in gray. The effect of resummation is roughly a 
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Table 3: Results for the total cross section in pb, using the default choice /i/ = 400 GeV. 
The first set of errors refers to perturbative uncertainties associated with scale variations, the 
second to PDF uncertainties. The most advanced prediction is the NLO+NNLL expansion 
highlighted in gray. 
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Table 4: Same as Table 3, but with the "educated" scale choice /// = rrit. 



10-15% enhancement over the fixed-order NLO result. A more important effect is that the 
resummation stabilizes the scale dependence significantly. Concerning the approximate NNLO 
schemes, the results from scheme A are noticeably higher than those from scheme B, but these 
differences are well inside the quoted errors. Since the two schemes differ only by terms 
proportional to 6{1 — z), this gives an indication of the size of the unknown constant terms. 

To some extent, the enhancement effect resulting from the resummation of the leading 
threshold terms can be mimicked using fixed-order results evaluated at a significantly lower 
factorization and renormalization scale /i/. Such an "educated" scale choice, which is often 
adopted in the literature on fixed-order calculations, is fif = rrit. Table 4 shows the cross- 
section predictions obtained in this case. The fixed-order results are indeed significantly 
enhanced with this scale choice. The resummed predictions, on the other hand, do not change 
much compared to those shown in Table 3. 
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Figure 13: Phase space in the (s, M^) plane. In the blue region along the diagonal threshold 
singularities arise, and the cross sections receives its main contributions. In the small green 
region near the origin Coulomb singularities appear and the small-/? expansion applies. 

6.4 Total cross section: Comparison with previous calculations 
6.4.1 SmalI-/3 expansion 

The approach pursued here offers an alternative to the direct threshold expansion of the total 
partonic cross section in the limit /3 — )■ 0, corresponding to s — )■ 4mf. In this case not only 
the phase-space for real gluon emissions shrinks to zero, but in addition the top and anti- 
top quarks are produced at rest in the partonic center-of-mass frame, which implies that in 
addition to soft-gluon singularities one encounters Coulomb singularities. The leading terms 
in the /3 — limit were first calculated at NNLO in [40,42], and later corrected in [45]. 

It is important to emphasize that the leading singular contributions to the total cross 
section arising from the /3 — ?■ limit do not coincide with those arising from the limit 2 — )■ 1, 
even after integrating over all kinematic variables. The reason is that after convolution with 
the PDFs there are no truly small scale ratios left in the process (the total center-of-mass 
energy ^/s at the Tevatron or LHC are so large that they can be taken to infinity compared 
with the scale m^). The large perturbative corrections to the cross section arise dynamically, 
because of the relatively strong fall-off of the parton luminosities combined with the fact that 
the partonic cross sections receive their dominant contributions from the region near Born- 
level kinematics [64,66]. One would then expect that the most accurate account of enhanced 
perturbative corrections should be the one that captures enhanced contributions in all relevant 
regions of phase space. 

In terms of the variables s and M^, the phase space is given by the triangular region 
Amf < < s < s, as illustrated in Figure 13. The large threshold terms considered in 
the present paper and in our previous work [56] are located along the diagonal, where the 
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partonic cross sections are largest. The large corrections arising in the /3 — )■ limit are located 
near the origin of the diagram, where both s and approach 4mf. The parton luminosities 
are largest for small values of s. It is obvious from the figure, and also by considering the 
invariant mass distributions shown in Figure 8, that the region near the origin gives only a 
very small contribution to the total cross section. This fact is most clearly demonstrated by 
the distributions in the variable Pt shown in Figure 11, which peak at Pt ~ 0.6 (Tevatron) and 
Pt ~ 0.7 (LHC). The region of small velocity, say below /3t = 0.2, obviously yields very small 
contributions to the total cross sections. Since the variable /3 is always larger than Pt, this 
conclusion is even more true for the small-/? region. On the contrary, the approach pursued 
in the present work accounts for enhanced perturbative contributions in all regions of phase 
space giving rise to large contributions to the total cross section. It is completely analogous 
to threshold (or soft-gluon) resummation for Drell-Yan or Higgs production at fixed value of 
the lepton pair or Higgs boson mass. Even though one can never be sure how accurately the 
full NNLO correction to a cross sections is approximated by a subset of calculable terms, we 
strongly believe that our treatment provides an approximation that captures more physics 
than that based on the /3 — )■ limit. We therefore expect our predictions for the invariant 
mass distribution, total cross section, and forward-backward asymmetry to be the most reliable 
available at present. 

Having just argued that the /3 — limit is not of much relevance for the total cross section, 
it is nevertheless interesting to study how well our predictions fare in this region. Since in our 
case the top quarks are generically not at rest in their center-of-mass frame, we are not dealing 
with Coulomb singularities, and hence our approximate prediction for the NNLO corrections 
to the cross section misses a subset of terms involving potential-gluon exchange. We will now 
study in more detail which of the singular terms in the /3 — )■ limit can be recovered in our 
approach. To this end, we write the total cross section in the form 

-;|E£f#«(;.'')/«(^.'')- (10^) 

We can obtain an expression for the perturbative functions fij by integrating (5) over all of 
phase-space, in which case we find 



We can now evaluate the above formula in the limit s — 4m^. Defining expansion coefficients 
for the functions fij as 



f^, = f\f+A^laJlf + {A^:as 



+ ..., (109) 



the answer for the scale- independent pieces with nh = 1, rii = 5, and = 3 can be written as 

910.22 In^ /3 - 1315.5 In^ /3 + 592.29 In^ /3 + 452.52 In /5 
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Table 5: Results for the total cross section in pb, using the default choice /U/ = 400 GeV. Some 
numbers from Table 3 are compared with results obtained from different implementations of 
the small-/? expansion (see text for explanation). The errors have the same meaning as before. 
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where the dots refer to (9(/31n°/3) terms, which are yet unknown. We have split the answer 
into the piece recovered from the expansion of our results, which we have written explicitly, 
and a piece related to NNLO effects from potential gluons, which would be recovered from 
the small-/? expansion of the as yet unknown /i-independent part of the NNLO coefficient 
Cq in (99). Such potential-gluon contributions were obtained in [45] by using the two- loop 
calculations of [89-91], and lead to the additional terms 

.potential ^ 3W7 ^ 1 ^504451^^ _ Q^27A) + 76.033 lu /3 , 

/5^ P 



^potential _ __ ^ _ (_79.270 lu /3 + 227.59) - 290.76 In /3 . 



The coefficients f-'^'^^ and f-f'^' (with ij = qq, gg) in (109) are also recovered from the expansion 



/32 /3 



:iii^ 



.(2,2) 



of our results, up to terms of 0(/3^). We have checked that they agree with the results given 
in [42] when expanded to that order. 

It is worth noting that after obtaining the small-/? expansion as in (110), one can replace 
the approximated Born prefactors, 7r/3/9 in the qq case and 77r/3/192 in the gg case, with 
the exact Born-level results. This procedure has been adopted in the recent literature on the 
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small-/? expansion [40,42,45]. Therefore, we will differentiate the two versions of the small- /3 
expansion: version 1 refers to the version where a complete expansion is carried out as in 
(110), while version 2 refers to the version with the exact Born-level results as prefactors. 

We are now ready to compare our results to those obtained using the small- /3 expansion. 
We focus first on a comparison at NLO, where the exact answers are known. In the upper 
portion of Table 5, we show the results of the different approximations as well as the exact 
NLO cross sections. As we have mentioned in the last section, our NLO leading singular terms 
are always a good approximation to the exact NLO results, at both the Tevatron and the LHC. 
On the other hand, the small-/? expansion version 1 tends to overestimate the cross section at 
the Tevatron by more than 2 pb. By incorporating the exact Born prefactors (version 2) the 
small-/? expansion works better, but still it overestimates the exact results by about Ipb. At 
the LHC, the small-/? expansion happens to give results closer to the exact answers. However, 
as we will now explain, this is a coincidence. In Figure 14 we plot the NLO corrections to the 
cross sections, including the parton luminosities, as functions of /?. The black solid curves show 
the exact results, the red solid curves our leading singular terms, and the dotted curves the 
results obtained using version 2 of the small-/? expansion (the dashed curves will be explained 
in the next subsection). In the small-/? region, all the approximations work rather well as 
expected. With increasing /?, the different approximations start to deviate from one another. 
We observe that, at both the Tevatron and the LHC, our approximations always reproduce the 
shapes of the exact results quite well, which is not at all achieved by the small-/? expansion. 
The fact that the small-/? expansion overestimates the cross section at the Tevatron, where the 
qq channel dominates, is evident from the left plots in Figure 14. At the LHC, where the gg 
channel dominates, the small-/3 expansion does not reproduce the shapes of the exact results, 
even though it happens that the integrated cross sections are close to the exact ones due to 
a coincidental cancellation. However, it is unlikely that a similar cancellation will happen at 
NNLO. 

We next compare our best prediction, NLO-I-NNLL, to the best prediction obtained using 
the small-/? expansion at NNLO, which is obtained by adding the approximate NNLO correc- 
tions derived using the small-/? expansion to the exact NLO cross sections. Without knowing 
the exact expression for the NNLO corrections, it is hard to tell which one is closer to the 
true answer, but we can study the validity of small-/? expansion by investigating the effects 
of the subleading terms in /? that are contained in our results. We have included in Table 5 
the numerical results for the cross section obtained by evaluating the small-/? expansion (110) 
of our approximate NNLO formula, without including the extra potential terms (labeled "/?- 
exp." in the table), and that obtained using this expansion plus the potential terms in (111) 
(labeled "/?-exp.-|-potential"). We notice that the NLO+NNLL resummed results and the 
small-/? expansion differ by about 10-15% at the Tevatron and about 6% at the LHC. On the 
other hand, the effect of adding the potential-gluon contributions to the small-/? expansion, 
which cannot be reproduced in our approach, is always smaller than 2%. We conclude that the 
bulk of the terms that become singular in the /? — )■ limit are accounted for in our approach. 
The reverse statement is not true. A resummation based on the /? — > expansion does not 
account for the bulk of the terms that become singular in the z ^ 1 limit, and our analysis 
suggests that subleading terms in /? are by no means generically small. 

The CDF collaboration has recently published a combined value for the ti cross section. 
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Figure 14: The corrections to da/dp for the different approximations mentioned in the text 
at the Tevatron and LHC, with /i/ = rrit. The plots on the left side show the qq channel, those 
on the right the gg channel. 

finding (Xcxp = (7.50 ± 0.48) pb [92]. They have compared their result with the theoretical 
prediction ath = (7.46!']ogQ) pb derived using the small- /3 expansion approach from [42]. This 
number is significantly higher than our best prediction ctnlo+nnll = (6-30;!;q;}91q;23) pb given in 
Table 3, but it is also higher than the value ctnnlo, /3-exp.+potentiai v2 = (6-95loi9la26) shown 
in Table 5, which is our own result obtained using the small-/? expansion. This latter difference 
is partly due to a smaller top-quark mass (172.5 GeV instead of 173.1 GeV), a different PDF 
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set (CTEQ6.6 instead of MSTW2008NNLO), and a different value of a,{Mz) (0.118 instead 
of 0.117), and partly due to a lower factorization scale {rrit instead of 400 GeV). Changing 
our input values to those used in [42], we would obtain o"nlo+nnll = (6.86lo;2olo;35) pb and 
cnnlo, /3-cxp.+potcntiai v2 = (^-SSlQi^ylasg) pb. We obscrvc that with these settings the predic- 
tions are indeed higher,^ and our (Tnnlo, /3-cxp.+potcntiai v2 is now very close to the result from [42]. 
The remaining slight difference can be attributed to the use of an incorrect two-loop coefficient 
in [42], which was later corrected in [45]. Our best prediction unlo+nnll is still about 10% 
lower than the result obtained using the small- /3 expansion, which can be understood from the 
discussions in the previous paragraphs, but is consistent with the experimental result within 
errors. 

6.4.2 Other approaches and kinematics 

Since in this paper we are interested in the invariant mass distribution of the tt pair, we have 
adopted the so-called pair-invariant mass (PIM) kinematics, which is defined by 

NliPl) + N2{P2) ti{ps + P4) + X{k) . (112) 

When calculating the transverse-momentum and rapidity distributions of the top quark, with 
the momentum of the anti-top quark integrated over, one instead considers the one-particle 
inclusive (IPI) kinematics, which can be written as 

iVi(Pi) + N^{P^) ^ t(p3) + X'[^{V,) , (113) 

where here the final state X' contains the t quark accompanied by additional emissions (see 
e.g. [50,51]). In both cases, the threshold limit corresponds to the limit in which these extra 
emissions are soft, implying — and p'l fnf, respectively. In the PIM case, this implies 
— )■ s, where = (^3-1-^4)^ and s = {pi + p2)'^ can be calculated without reference to 
the unobserved momentum k. In the IPI case, it implies 84/3 = (1 — 2E^/\/§) 0, where 
Si = s + ti + ui can be calculated from (3) without reference to the unobserved momentum 
P4, and -E3 denotes the energy of the top quark in the partonic center-of-mass frame. 

Although the PIM and IPI kinematics are applicable in different differential distributions, 
they can both be integrated over to obtain the total cross section. In the sense that they 
are both applicable in the limit where the extra emissions are soft, threshold resummation 
for the total cross section based on PIM or IPI kinematics amounts to resumming the same 
leading contributions, but differs by subleading corrections. Several authors have found signif- 
icant numerical differences between the results obtained using the two kinematic schemes (see 
e.g. [50-52]). Therefore, it is interesting to work out the IPI kinematics also in our approach. 

Before moving onto IPI kinematics, it is however necessary to point out an important 
difference between our results and previous ones obtained using PIM kinematics. In the tradi- 
tional approach, the leading singular terms in {1 — z) are written in terms of the distributions 

^The difference between the central values of our NLO+NNLL predictions (6.86 pb vs. 6.30 pb) is 0.56 pb, 
of which about 0.18 pb is due to the scale choice, about 0.10 pb is due to the smaller mt, and another 0.28 pb 
is due to the different PDF set and as{Mz) value. 
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P„(^) = [ln"(l-z)/(l 
of the distributions 



z)]+, while in our approach they are more naturally written in terms 



I- Z 



M^(l-z) 



(114) 



The additional factor of z in the logarithms is a subleading effect, but it is relevant in practice. 
This has been studied in detail for the simpler cases of Drell-Yan [64] and Higgs production [66], 
but analogous remarks hold also in the present case. Of crucial importance in this context is the 
fact that our resummation method works directly in momentum space [62], and as a result the 
matching onto analytical fixed-order expressions is particularly transparent. We have observed 
in previous work [64, 66] that at NLO in perturbation theory all logarithmic corrections to 
Drell-Yan and Higgs production, even those entering at subleading orders in (1 — z), can be 



expressed in terms of the "hard" and "soft" logarithms ln(M^/ /i^) and ln(s 

\2 



-zY / fi'j), where 

= M^+i_ or m'jj as appropriate, and s = M'^/z. Choosing s (1 — z)^ = M^(l — z^/z as the 
argument of the soft functions in the SCET approach thus provides an "optimal" matching 
onto fixed-order expressions, which automatically accounts for subleading logarithms In z. We 
have observed that, with this choice, the leading singular terms near partonic threshold account 
for the bulk of the NLO and NNLO corrections to the Drell-Yan and Higgs production cross 
sections. In the present paper, we have shown that the same is true for the NLO corrections 
to the ti production cross sections at the Tevatron and LHC Our resummation procedure is 
thus different, in the way we treat subleading terms, from the traditional PIM scheme, and 
for the reason just explained it is bound to be more accurate. This can be clearly seen from 
Figure 14, where we have plotted the NLO leading singular terms both in our scheme (red 
solid curves) and in the traditional PIM scheme (dashed curves). One can find that our results 
always provide a better approximation to the exact NLO correction than the traditional PIM 
results. Especially for the gg channel at the LHC, where the traditional PIM results were 
argued to be a bad approximation [51,52], our result still provides a decent approximation. 
We now turn to the IPI kinematics, in which case the relevant differential cross section is 



dTi dUi 3s2 



E 



dxi dx2 

Xi X2 



fi/N^ (Xl, /i/) fj/N:,{x2,fif) C[j{Si, ti, Ui, mt,flf) 



;ii5) 



Here Ti = {Pi-psf-mf and Ui = (P2-P3 



-m 



are the hadronic counterparts of ti and ui. 



Near the partonic threshold S4 — > 0, the hard-scattering kernels C^' • can again be factorized as 



;il6) 



where the hard functions Hij are the same as in the PIM case, once we re-express and 
cos^ in terms of ti and Ui. The new soft functions S^j can be calculated from (68) with the 
replacement — ?■ ■ k. After these changes, all techniques we have developed in this paper 
can be applied to the IPI case, too. In particular, we can derive a resummed expression for 
the hard-scattering kernels 



C'{s4,ti,ui,mt,fif) = exp 



2ar(/is,/i/)ln— ^ + 4a^0(/i^, /i/) 

tiUi 
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X Tr 



U{ti, ui, mt, Hh, l^s) H{ti, ui, mt, Hh) U\ti, ui, rrit, Hh, l^s) 



S4 



S4 



r(2r/) 



;il7) 



By setting fih = fis = fJ^f the above formula, we can derive the leading singular terms in 
IPI kinematics. We have checked that the leading terms at NLO agree with those in the 
literature [14,59]. A numerical analysis of the resummation in IPI kinematics is better to 
be done in the context of transverse-momentum and rapidity distributions of the top quark, 
where a systematic study as we have done for PIM kinematics can be carried out to determine 
the optimal choices of the scales fih, fig and /i/. Therefore we will leave it for a future analysis. 
Nevertheless, given the arguments in the previous paragraphs, we do not expect significant 
differences in the results for the total cross sections obtained using the two kinematic schemes. 



6.5 Charge and forward- backward asymmetries 

The top-quark pair charge asymmetry is an important observable at the Tevatron, which 
originates from the difference in the production rates for top and anti-top quarks at fixed 
scattering angle or rapidity [20,21]. To describe the total charge asymmetry, we first introduce 
a charge-asymmetric cross section as the difference between the production of top and anti-top 
quarks in the forward direction. 



Aar 



1 rV^ 
dcose / dM 

J2mt 



dMd COS e dMd cos 9 



:il8) 



The total charge asymmetry is then defined by A^. = Aac/cr, where a is the total cross 
section. As a consequence of the charge conjugation invariance of QCD, Aac is the same as 
the difference between the cross sections for the production of top quarks in the forward and 
backward directions, 



Aar = Aa 



FB 



dM 



2mt 



rfcost' 







dMd COS 9 



a cosy 

1 



dMd cos 9 



(119) 



Therefore, the charge asymmetry can also be interpreted as a forward-backward asymmetry 
for the top quark, defined as Ap^ = Actfb/c"- 

At leading order in QCD the charge-asymmetric cross section Aac vanishes. This 

quantity receives non-zero contributions starting at 0{a^). These arise if, in the interference of 
one-loop and tree-level diagrams, the top-quark fermionic line and the light-quark fermionic 
line are connected by three gluons. The same observation applies also to the interference 
of two tree-level diagrams with three particles in the final state. In Figure 15 we show the 
interference of the planar box with the tree-level diagram and the corresponding interference 
of real emission diagrams. The other contribution to the asymmetry at O(a^) in the quark- 
annihilation channel originates from the interference of the crossed box and tree-level diagram 
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Figure 15: Examples of interferences contributing to the charge-asymmetric cross section. The 
two-particle cut corresponds to the interference of a one-loop box with the tree-level diagram, 
while the three-particle cut corresponds to the interference of tree-level diagrams with a ttg 
final state. 

(or from the corresponding real emission case). This can be visualized by imagining to cross 
the two gluons on the left side of the heavy-quark triangle in Figure 15. The color factors 
multiplying the structure in Figure 15 or its crossed counterpart are respectively 

Cplanar = Jg^y^ {fabc + '^abc) ' ^crossed = Jg^y^ i'^abc ~ fabcj ' (-'-^O) 

where /f^^ = (A^^ — 1)A^ and c?^^^ = (A^^ — 1)(A^^ — 4)/iV. When the color factors are stripped 
off, the interference in Figure 15 and its crossed counterpart satisfy the relation 

= -c^f^™? • (121) 

This relation holds both for the three-particle and the two-particle cuts. Therefore, the charge- 
asymmetric cross section is proportional to d\^^. The interference of the one-loop box diagrams 
with the tree-level diagram gives a positive contribution to the asymmetry, which is partially 
canceled by the asymmetry originating from the interference of initial- and final-state ra- 
diation diagrams. An additional small contribution to the asymmetry at 0{al) originates 
from the flavor excitation channel gq{q) — ?■ tiX at tree level. The gluon-fusion channel does 
not contribute to the charge-asymmetric cross section. The study of the charge-asymmetric 
cross section at 0{al) shows that top quarks (anti-top quarks) are preferably emitted of the 
direction of the incoming quark (anti-quark); consequently, in pp collisions top quarks are 
preferably emitted in the direction of the incoming proton [20,21]. 

In the asymmetric cross section (118), the definition of the scattering angle 9 depends on the 
reference frame, therefore the value of the asymmetry is also frame-dependent. Experimentally 
the forward-backward asymmetry has been measured both in the laboratory frame and in the 
tt rest frame [10, 11], where the result in the former frame is about 70% of the result in the 
latter. Near the partonic threshold the ti rest frame and the partonic center-of-mass frame 
coincide. Using the complete 0{al) expressions for Actfb in [21], we have checked explicitly 
that the results in these two frames only differ slightly. It is therefore interesting to employ 
the formulas derived in the previous sections to calculate the forward-backward asymmetry 
in the partonic center-of-mass frame. The resummation of partonic threshold effects in the 
forward-backward asymmetry at NLL order in Mellin space was first considered in [49]. In 
what follows we extend these calculations to NNLL order, and also study the effect of the 
approximate NNLO corrections in fixed order. 
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0.2 < /i//TeV < 0.8 


mt/2 < fif < 2mt 


AapB [pb] 


^*FB [%] 


Ao-pB [pb] 


^*FB [%] 


NLL 

NLO, leading 
NLO 


n iq+0.09 


5.811^2 

c 9+0.4 

« 7+0.6 
' -0.4 




U-'J-L-O.IO 


c q+3.4 

5.7l°i 
7.41°:^ 


NLO+NNLL 


n 40+0 06 




Q 4c;+0.08 


' •'-'-0.7 


NNLO, approx (scheme A) 
NNLO, approx (scheme B) 




a 4+0-9 

c; O+0.8 
^■0-0.6 


48+° " 

n 45+0.09 


7 5 + 1.3 
' •'J-0.9 

^•'-'-0.8 



Table 6: The asymmetric cross section and forward-backward asymmetry at the Tevatron, 
evaluated at different orders in perturbation theory in the partonic center-of-mass frame. The 
errors refer to perturbative uncertainties related to scale variations, as explained in the text. 

The results of our analysis are summarized in Table 6. In the column labeled Actfb; we 
show the total charge-asymmetric cross section evaluated at different orders in perturbation 
theory. As in the previous sections, we have used rrit = 173.1 GeV and adapted the order of 
the PDFs as described in Table 2. For the factorization scale, we consider our standard choice 
/i/ = 400 GeV as well as fif = rrit. The results in resummed perturbation theory use fih = M 
and fig chosen according to (105) by default, and the uncertainties are obtained by varying 
these up and down by a factor of two and adding the different uncertainties in quadrature. 
The uncertainties in the fixed-order results are obtained by varying /.i/ up and down by a 
factor of two. The counting used in the table refers to the order at which the differential cross 
section itself is needed, relative to a^; this differs slightly from the counting in fixed order 
used in [21], which would count, for instance, our NLO as LO. To obtain the result in fixed 
order at NLO, we have used the formulas in Appendix A of [21]. At both /i/ = 400 GeV and 
/If = rrit, the NLO threshold terms recover about 80% of the full result in fixed order, which 
is roughly in line with our findings for the NLO corrections to the cross section. The table 
also includes our results for the forward-backward asymmetry. In calculating the asymmetry, 
we first evaluate the numerator and denominator of the ratio Ap^ = AapB /o" to a given order 
in RG-improved or fixed-order perturbation theory, and then further expand the ratio itself. 
When performing the calculation in this way, the errors in the asymmetry at NLO are actually 
smaller than those at NLO-I-NNLL order, even though the scale variations in the numerator 
and denominator of the ratio are much larger. We note, however, that if we chose instead to 
not further expand the ratio, the NLO result would be decreased by about 20%, while the 
NLO-I-NNLL result changes only by about 5% and should therefore be considered the more 
reliable prediction. In [21], an overall factor of 1.09 due to electroweak corrections is included 
for the asymmetry. However, a smaller correction was recently obtained in [93], and it is also 
scale-dependent. Therefore we have chosen not to include these corrections in our results. 

Our results are in good agreement with the previous findings reported in [49,94]. We 
have not performed the calculation in the pp frame, but expect that the boost to this frame 
decreases the asymmetry by roughly 30%, as found in [94]. On the other hand, the CDF 
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Figure 16: Dependence of the total cross section on the top-quark mass defined in the pole 
scheme. The NLL and NLO+NNLL bands indicate the scale uncertainties and PDF uncertain- 
ties combined linearly, while the blue band represents the dependence of the DO measurement 
of the cross section on rrit, as described in [100]. For comparison, a CDF combined measure- 
ment [92] is also shown. 

collaboration at the Tevatron recently reported the value Ap^{exp) = (19.3 ± 6.9) % for the 
asymmetry in the pp frame [12]. The measured asymmetry exceeds the predicted one by about 
two standard deviations. In light of our results, we conclude that higher-order QCD effects 
are not sufficient to explain the large experimental value. Possible explanations of this fact in 
the framework of several new physics scenarios were recently investigated in [94-99]. 

6.6 Sensitivity to the top-quark mass 

Up to this point our results were obtained using the default value mt = 173.1 GeV for the 
top-quark mass defined in the pole scheme. Figure 16 shows the dependence of the total cross 
section on the value of m^. We show our NLL and NLO-I-NNLL results as bands representing 
the linearly combined errors from the scale uncertainties and PDF uncertainties, which has 
been converted to 68% CL (Icr) to match the confidence level of the experimental errors. We 
compare our result with a recent combined measurement of the total cross section performed by 
the DO collaboration [100]. We also show a combined measurement by the CDF collaboration 
[92]; however, since they did not provide the mass dependence of the cross section, we cannot 
use their results to extract the top-quark mass. From the overlap region of the bands we 
determine the top-quark mass to be mt = {163.0lil'l) GeV, which is lower than the latest 
world average nit = (173. l]*^]^ !^) GeV [1]. One should, however, interpret this "discrepancy" 
with care. It is unclear if the quoted world-average top-quark mass really refers to the pole 
scheme, and even if it does, to which order of perturbation theory it corresponds. Since the 
pole mass of a quark is not well-defined beyond perturbation theory, the value of rrit in the pole 
scheme changes significantly in different orders. It would therefore be preferable to perform 
both the theoretical and experimental analyses with a more physical mass parameter such as 
that defined in the MS scheme (see below). 
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Figure 17: Position of the peak of the invariant mass distribution (left) and of the mean 
invariant mass (M) (right) as functions of the top-quark mass. 





Mpeak [GeV] 


(M) [GeV] 


LO 


375.6l°i 




NIL 


374.5l°;2 




NNLL 


376.lt°i 


445.5^°:^ 


NLO+NNLL 


375.2l°;3 


445.0^°:^ 



Table 7: Values of the peak position of the mean invariant mass, obtained with rrit = 173.1GeV. 



As we have seen, there are still rather large uncertainties in using the total cross section as 
a mean to extract the top-quark mass. We therefore investigate what additional information 
can be gained from the invariant mass distribution. On the left in Figure 17, we show the 
position of the peak of the invariant mass distribution, Mpeak, as a function of rrit. On the 
right we show the corresponding dependence of the mean invariant mass (M) . We do not show 
the errors arising from scale variations, as they would be too small to be visible on the scales 
of the plots. As illustrated in Table 7, both observables turn out to be very stable under the 
change of the order of perturbation theory. (Observing that the shifts between different orders 
are larger than those indicated by the scale variations, we however note that scale variation 
might not be a good estimator of the uncertainty.) They show an almost linear dependence 
on the mass of the top quark, which leads to a nearly linear translation of errors. A precise 
measurement of one of these observables would thus lead to a direct measurement of rrit with 
about the same precision. Of course this can only be taken as a qualitative remark at the 
moment, but with increasing data this might become an important strategy for a precision 
measurement of rrit. 

It is a well-known fact that the pole mass defined in perturbation theory is an ill-defined 
concept, as it is plagued by renormalon ambiguities in higher orders of perturbation theory 
[101, 102]. While the corresponding intrinsic uncertainty in rrit is of order Aqcd and thus 
smaller than the present experimental error in this parameter, it is nevertheless useful to 
consider predictions for the cross section (and distributions) parameterized in terms of a more 
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mt [GeV] 


LO 


164.0 


NLO, NLL 


171.7 


NNLO, NNLL 


173.1 



Table 8: Values of the pole mass rrit in different orders of perturbation theory, corresponding 
to a fixed value rnt{mt) = 164.0 GeV. 



suitable mass parameter. For the purpose of illustration, we investigate in the following the 
impact of using the MS mass fnt{fi). We do this by changing the value of the pole mass 
at different orders in perturbation theory, such that the "physical" mass mf(/i) remains the 
same. This is implemented through the relation 



(122) 



where for QCD (A^ = 3) with Uf = 5 active, massless flavors [103-105] 



,m 4 ^ 2053 tt^ tt^ , (3 379 , 37 . , , 

= - + L„ , = + — + — n 2 - — + L^ + —LL, 123 

3 ' 288 18 9 6 72 24 ^ ^ 

with Lm = ^^{fJ'^ /^tif^))- ^se a flxed input value frit = mtijnt) = 164.0 GeV, chosen such 
that at two- loop order the corresponding pole-scheme parameter mt in (122) coincides with 
our default value mt = 173.1 GeV. In lower orders we adjust mt such that the MS mass stays 
unchanged, which leads to the values collected in Table 8. In Figure 18, we show the impact 
on the invariant mass distribution of choosing mt in the way described above. Compared with 
Figure 8, we observe an improved convergence of the perturbation theory, both in flxed order 
and after threshold resummation. This flnding, which has previously been made in [42], could 
have been anticipated by looking at Figure 16 and Table 8. The increase of the cross section 
and invariant mass distributions in higher orders is to a large extent compensated by the effect 
of increasing the value of the pole mass. 



7 Conclusions 

We have studied higher-order perturbative corrections to the differential cross section for tt 
hadroproduction in the partonic threshold region, where the invariant mass of the tt pair 
approaches the partonic center-of-mass energy. This involved using SCET to derive a fac- 
torization formula expressing the hard-scattering kernels as products of matrix-valued hard 
and soft matching coefficients in this kinematic regime, calculating these coefficients to NLO 
in perturbation theory, and solving the RG equations for these functions directly in momen- 
tum space to achieve NNLL accuracy for the resummed differential cross section. We also 
presented an approximate NNLO formula, which includes all terms proportional to singular 
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Figure 18: Fixed order (left) and resummed (right) results for da/dM, changing rrit at different 
orders in perturbation theory as shown in Table 8. 



plus distributions in the variable {1 — z), as well the scale-dependent part of the coefficient of 
the (5(1 — z) term. The momentum-space resummation techniques allow for a straightforward 
matching of the resummed results with the exact results at NLO in fixed-order perturbation 
theory. We used this fact to perform numerical studies of the invariant mass distribution, the 
total cross section, and the forward backward asymmetry at NLO+NNLL order. 

Our results are expected to work progressively better when the invariant mass M be- 
comes large, since in that case the partonic threshold region becomes more important due to 
kinematics. However, our numerical results show that the full result at NLO in fixed-order 
perturbation theory is dominated by the leading terms in the threshold expansion, even at 
low values of invariant mass near the peak of the distribution. This provides evidence for a 
dynamical enhancement of the partonic threshold region. Assuming this to be the case also at 
higher orders in perturbation theory, we have argued that our results obtained at NLO-I-NNLL 
order provide the best available predictions for the ti invariant mass spectrum in the whole 
range of M. The benefits of performing the resummation are best seen through the reduction 
of scale dependence and the improvement of the convergence shown in Figure 7-10. Compared 
to the fixed-order NLO results, the RG resummation reduces the scale uncertainties in the 
differential cross sections and therefore provide the more reliable predictions. 

In view of the accuracy we can achieve for the differential cross section, we have also made 
predictions for the total cross section and the forward-backward asymmetry by integrating 
the differential distribution over the appropriate regions of phase-space. We believe that 
these results, summarized in Tables 3, 4 and 6, are currently the best predictions for these 
observables. In particular, we have compared our results for the total cross section with 
those obtained in the conventional approach of expanding around /3 — )■ 0, where the top- 
quark pair is produced nearly at rest. By expanding our approximate NNLO formulas in the 
limit /3 — > and dropping sub leading terms in /3, we recover the results recently obtained 
in [45], up to terms at NNLO related to potential-gluon exchange. A numerical analysis 
showed that the contributions of these extra potential-gluon corrections are actually smaller 
than the subleading terms in /3 contained in our formulas, whose contribution is about 6-15%. 
This is not surprising, because the cross section receives the dominant contributions from 
regions with /3 ~ 0.6, which is not especially small. For this reason, it seems to us that the 
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results obtained here offer a legitimate alternative for including higher-order corrections to the 
total cross section. As for the forward-backward asymmetry, our results are in quantitative 
agreement with previous findings, and thus cannot explain the large value measured by the 
CDF collaboration. 
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A RG-evolution factors and anomalous dimensions 

Here we collect the expressions of the perturbative solutions to NNLL order for the RG factors 
S, a^, and u in (83) and (84), which appear in the solution of the evolution matrix U in the final 
resummed result (92). We also collect the expansion coefficients of the anomalous dimensions 
and QCD /3-function needed at this order. 

We first define expansion coefficients of the anomalous dimensions and QCD /3-function as 



(3 (as) = -2as 



(A.l) 



and similarly for the other anomalous dimensions (recall that Fcusp = C*F7cusp for the qq 
channel, and Fcusp = C'aTcusp for the gg channel). In terms of these quantities, the function 
Or is given by [63, 106] 

and the result for the Sudakov factor 5* reads 
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where r = as{fi)/as{h'). 

The matrix equation for u can be evaluated using the techniques described in [107, 108] 
To this end, we first define the matrix V which diagonahzes 7'^(o) as 



(A.4) 



We also define the vector 7^^'^) consisting of the diagonal elements of 'y'l^^^ ■ Then the solution 
at NNLL reads 



u{M,mt,cose,fih,fi) = V + 



2/9o 



D 



An 



where the matrix elements of K are 
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Finally, we collect the expansion coefficients of the anomalous dimensions and QCD /3- 
function needed in this work. The cusp anomalous dimension to three-loop order is [109] 
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For the anomalous dimensions entering 7^^, we have [74,79] 
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7o^ = -2Cp , 

= CfCa {-- + — - ^C- 
as well as [54,55,74,110-112] 
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The anomalous dimensions of the PDFs are given by 
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and for the QCD /3-function to three-loop order we have 
/3i = yCl - yC^T^nj - ACfTfu^ , 
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